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ABSTRACT 

We present comprehensive multiwavelength observations of three gamma-ray bursts (GRBs) with 
durations of several thousand seconds. We demonstrate that these events are extragalactic transients; 
in particular we resolve the long-standing conundrum of the distance of GRB 101225A (the "Christmas- 
day burst"), finding it to have a redshift z = 0.847, and showing that two apparently similar events 
(GRB 111209A and GRB 121027A) lie at z = 0.677 and z = 1.773 respectively. The systems show 
extremely unusual X-ray and optical lightcurves, very different from classical GRBs, with long lasting 
highly variable X-ray emission and optical light curves that exhibit little correlation with the behaviour 
seen in the X-ray. Their host galaxies are faint, compact, and highly star forming dwarf galaxies, 
typical of "blue compact galaxies". We propose that these bursts are the prototypes of a hitherto 
largely unrecognized population of ultra-long GRBs, that while observationally difficult to detect may 
be astrophysically relatively common. The long durations may naturally be explained by the engine 
driven explosions of stars of much larger radii than normally considered for GRB progenitors which are 
thought to have compact Wolf-Rayet progenitor stars. However, we cannot unambiguously identify 
supernova signatures within their light curves or spectra. We also consider the alternative possibility 
that they arise from the tidal disruption of stars by supermassive black holes. 

Subject headings: galaxies: distances and redshifts - gamma rays: bursts - techniques: photometric 
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1. INTRODUCTION 

Classical long-duration gamma-ray bursts (GRBs) are 
active over timescale^^ ranging from ~ 2 s up to several 
hundred seconds in the observer frame. In the context of 
the conventional collapsar model, this period is thought 
to reflect the initial lifetime of the ultra-relativistic jet, 
which, once it has drilled through to the surface of the 
stripped-envelope progenitor star, goes on to p roduce the 



^^prompt" p hase of high energy emission (e.g., Bromberg 



etaI]|2QT2| . 

Over the years a small number of apparently 
even longer events have been observed, for exam- 
ple GRB 970315 had tgo ~ 1360 s, and was one of 
several bursts (out of ~ 3000) with tgo > 500 s 
discovered by the Compton Gamma-Ray Observa- 
tory (CGRO) /B mst and Transient Source E xperiment 
(BATSE) (e.g., [Tikhomirova fc Stern] |2QQ5| , although 
in several cases this very long duration was due to the 
tail of emission for extremely bright fast ri se exponentia l 
decay bursts, such as GRB 971208 (Giblin et al.'"2002"). 
Other instruments have also seen very long events, ex- 
amples i nclude the 1600 s GRB 020410 detected by Bep- 
poSAX (Nicastro et~ar]|2004| [Levan et ap '2005|,_and 
GRB060814B detected by"KQnus-Wind ([PaPShin et al.| 
|008). For ah missions, such bursts appear to be ex- 
tremely rare events, although the ability of a given in- 
strument to find such bursts is a strong function of the 
detector characteristics and trigger algorithms as well as 
the observing scheme adopted. Swift has only seen a 
handful of GRBs with com parably lo ng durations . In 
two cases, XRF 060218 (e.g., iP ian et al. 2006; Campana] 
|et aL'^^OOe"; "Soderber g etjd^ ^ and XRF 100316D 
(e.g.. Starlin g et al. 20TT| " (Jhornock et al.^^2010b| ), these 
turned out to be low-redshift, low- luminosity explosions 
whose nature appears distinc t from the bulk of c osmo- 



logical GRBs (e.g., [So derber ^ et al | |2004l |2006[ [Liang 
|et al."2007 Chapman et al.||2007[ ), although such events 



provide some of the best evidence for the as sociation of 



GRBs with broad-lined t ype Ic super novae (jPian et al. 
2QQ6||Bufano et al.||2012D . 



Interest in the nature of very long GRBs has been 
brought to the fore by a series of recent discoveries 
by Swift. GRB 101225A (the so called "Christmas-day 
burst") was active w hen Swift first slewed to the field 
( Palmer et al.||20l'Qa), and remained s o for several thou- 
sand seconds longer ( Thone et al. 2011 ). Indeed, its prop- 
erties were sufficiently unusual that two entirely different 
models for its origin have been proposed. In the first, the 
burst arises from the tidal shredding of an asteroid by a 



neutron star within our own Galaxy (Campana et al. 



2011), while the second posits an extragalactic burst. 



associated with a supernova inside a dense envelope at 
z = 0.33 dThone et ar||201lD . The obstacle to further 
progress was the lack of a robust determination of its 
distance. 

In another case, that of GRB 110328 A/ Swift 
J1644+57, the prompt phase spanned several days 
and has been argued to be most likely a tidal disruption 
event (TDE) produced when a star is consumed by a 

■^^ Conventionally GRB duration is most often quantified as the 
period, tgn, over which 90% of the gamma-ray emission is observed 
(e.g.,[Kouveliotou et al. 1993). Strictly speaking this is dependent 
on the sensitivity and bandpass of the detecting instrument. 



et al. 2011 Burrows et al 



supermassive black hole (e.g. [ Levan et alT 2011| Bloom 



Zauderer et al.||2011b). 



An important piece of supp orting evidence tor this 
interpretation was the precise spatial coincidence of the 
outburst with the nucleus of its host galaxy (Levan 

2011b). However, it should 



et al. 



2011 



Zauderer et al. 



be noted that this event has provoked lively debate 
and some collapsar-like mo dels have been proposed 



as alternative explanatio ns ([Quataert fc Kasen 2012 
[Woosley fc Hegerl 2012); the total energetics being 



comparable in both scenarios (i.e., around the rest-mass 
energy of a star). 

In this paper we present evidence for another distinct 
class of ultra- long duration GRBs. We analyse new ob- 
servations and spectroscopy of GRB 101225A, the proto- 
type of this class and use extensive multiwavelength ob- 
servations to show that GRB 11 1209 A is similar (see also 
the recent paper by Gendre et al. 2012). Other members 
of the same class may also be present within the Swift 
sample, and we discuss the recent GRB 121027A as a 
likely example of such an event. The paper is structured 
as follows: we investigate the high-energy properties and 
X-ray lightcurves and spectra of the bursts in section 
[2j and then present spectroscopy and imaging to estab- 
lish distances (section [3|, optical lightcurves (section [4| 
and host galaxy properties (section p|; from this we con- 
sider constraints that can be placed on their progenitors. 
Throughout we focus particularly on GRB 101225A and 
11 1209 A, since these are at a lower redshift and are bet- 
ter studied (section [t]) . Finally in section [s] we consider 
if other similar examples might exist within the archival 
Swift population, but have been unrecognized to date. 

2. HIGH ENERGY PROPERTIES 
2.1. Prompt emission 

As discussed above, when initially discovered, both 
GRB101225A and GRB111209A stood out as being of 
remarkably long duration. In fact, in each case the detec- 
tion by Swift was made as an "image trigger" ; in other 
words the flux detected by the Burst Alert Telescope 
(BAT) did not itself rise above the rate threshold^ but 
instead a statistically significant source was founam the 
reconstructed image plana^ It is worth noting that in 
the case of GRB 11 1209 At his image trigger criterion 
was actually reached on two separate occasions, making 
it one of very few bursts t hat have triggered the instru- 
ment multiple times (Palmer et al.pOllp^ 

In the case of GRB 121027A, the prompt emission was 
sufiicient to provide a st andard "rate" trigger of the BAT 
(Barthelmy et al. 2012), but continued at a lower level 
for several thousand seconds beyond this point (Starling 
et al. in prep). 

One consequence of the differing modes of detection 
is that the trigger-time becomes particularly poorly con- 
strained for image triggers. The instrument only regis- 
ters a detection at the end of the image trigger period. 

Formally, there is not a single rate trigger, but several different 
rate-based triggering criteria. 

BAT is a very wide field coded-mask instrument. Most GRBs 
are initially identified by a rise in the overall count rate in the 
instrument ( jGehrels et al.|2 004 ). 

29 The very long event GRB 110S28A/ Swift J1644+57 also trig- 
gered the BAT on several occasions in this case over the space of 
~ 2 days 
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which can be significantly after the onset of bursting ac- 
tivity. Determination of the early power-law decay rate 
of the light-curve, for example, depends on the assumed 
start-time of the event, and so will be less certain for 
these events than more typical GRBs. For all bursts con- 
sidered in this paper we utilise trigger-times as reported 
by Swift, but note that these have significant uncertain- 
ties associated with them. 

Estimating the duration of very long events is sim- 
ilarly not straightforward. Firstly, they are active for 
much longer than the Swift orbit of ~ 90 minutes (and 
at fiux levels too faint for many other missions), and 
so there are periods of activity that are simply missed. 
Since GRBs are intrinsically highly variable this means 
it is rarely possible to reconstruct the missing fiuence 
over these gaps. This is important in these cases since 
the fiuence is dominated not by the bright peaks at the 
beginning of the lightcurve (indeed ultra-long bursts fre- 
quently do not have such peaks), but by the integration 
of the lower-luminosity, but much longer lived emission. 
In practice this means that, especially in the absence of 
data from other observatories that do not suffer from 
Earth occultation, any attempt to determine a tgo will 
necessarily be crude. It should be noted that Swift has 
recently introduced a new trigger based on integrated 
fiuence, which may increase the observed rates of these 
events 

In the case of ORB 101225A, we face the additional 
problem that the source was active when Swift first 
slewed to the location, so at best we have a lower limit 
on the duration. It was also clearly ac tive in the sub- 
sequent Swift orbit (Thone et al. 2011). This suggests 
a duration in excess ot -7000 s. l^br GRB 111209A 
there was fortunately also a detec tion from the Ko nus- 
WIND satellite in waiting mode ( Gole netskii et al.|2011 ) . 
This lightcurve clearly shows significant structure over 
a period of — 10000 s, including fiaring activity while 
Swift was in Earth occultation. Finally, in the case of 
GRB 121027A, the burst is clearly detected in subsequent 
Swift orbits, at a rate a factor of 5 lower than the early 
peak (Starling et al. in prep.), again suggestive of a dura- 
tion in excess of —6000 s, while observations with MA XI 



emission CTho ne et aLlMTj [Golenetskii et al.|2011|[G^ 



show ongoing activity during the Swift orbit gap (Serino 
et al.|[2Q12 ). We adopt the aforementioned durations as 



indicative of tgo, but note that they suffer from signifi- 
cant uncertainty. 

The approximate location of these events in the 
spectral-hardness versus duration plane is shown in 
Fig. [l} and in the average-luminosity versus duration 
plane in Fig |2j Both views provide good illustrations 
of their extreme natures with respect to other GRB and 
GRB-like populations. 

2.2. X-ray Lightcurves 

The X-ray (luminosity versus rest frame time) light 
curves of all three GRBs are shown in Figure [3j All 
exhibit long-lived high luminosity activity, punctuated 
by large scale dipping and fiaring, followed by a rapid 
decline roughly — 10^ s from the trigger. The bursts 
are detected in 7— rays during this plateau, and appear 
to show a broad correlation between the 7— and X-ray 

www. astro . Ijmu. ac .uk/grb2 012/presentations/| 

|presentations/g ehrels_liverpool2012 . pdf 



dre et al. 2012), albeit that the 7-ray detections are at a 
lower signal-to-noise ratio. In the case of GRB 101225A 
and GRB 11 1209 A the lightcurves are strikingly similar, 
while GRB 121027A exhibits rather higher magnitude 
variability. The observed behaviour is conspicuously un- 
like that seen in most GRBs (also shown), which typi- 
cally fade rapidly on timescales of hundreds of seconds, 
followed by more slowly fading afterglow emission. While 
late time activity is not especially unusual in GRB after- 
glows, manifesting as either fiares (normally within the 
first t housand seconds) or long-lived plateaux (e.g.jZhan^ 



et aLjflOOGj), the longevity of the emission at particularly 



high flux levels does make these bursts stand apart, as 
clearly seen in Figure [3) Below, we briefiy discuss the 
quantitive properties ofthe ultra-long burst light curves. 

For GRB 101225A, the lightcurve can be fit by a simple 
broken power-law (of the form F{t) ex t~"), with pre- and 
post-break indices of ai = and a2 = 5.86^0 195 

and a break at tb = 22000 ± 300 s, although the pro- 
nounced dipping behaviour means that the quality of the 
fit is relatively poor (x^/dof ~ 3.6) and, given this, the 
error bars reported on fit parameters should be treated 
with caution. We note that there is a hint of possible 
periodicity in the GRB 101225A light curve. Specifically, 
there is a sharp rise at the beginning of the second orbit 
of observation, and a dip at the end of it. Similarly, there 
appear to be dips just at the end of the third and fourth 
orbits. This could be indicative repetitive behaviour with 
a period around the Swift orbital period (or half of it), 
however, given the small number of dipping features and 
the lack of any apparent periodicity in the light curves 
(7-ray or X-ray) of the other two ultra- long bursts, we 
consider it most likely to be just coincidental in this case. 

For GRB111209A, the fits prefer a power-law with 
three breaks, although again the fit statistics are signifi- 
cantly impacted by dipping and fiaring behaviour (with 

X^/dof =7.5 in the final fit). In this case, we find the 

+0.01 

-0.01' 



relevant slopes to be ai = 0.28^0.015 ^2 = 0.77; 



as 5.30 



+0.06 
-0.06 



and 0^4 = 1.361qq5, ^^^^ break times 
of tbi = 2050 ± 20 s, tb2 = 16300 ± 100 s and tbs = 
46000 ± 1000 s. Although differing in details, the broad 
properties are very similar between the two bursts, in 
particular the timing, and presence of the steep decay 
index at times > 10^ s and with slopes steeper than 
a - 5. We note that in the case of GRB 11 1209 A the X- 
ray lightcurve returns to a typical afterglow decay value 
at the end of the steep slope. No such component is 
seen for GRB 101225A, although this may be due to the 
observations being of insufficient depth. 

The lightcurve of GRB 121027A is somewhat different: 
it appears to show an early fiare with a rapid decline fol- 
lowed by a very rapid renewing of the activity around 
1000 s after the initial trigger, persisting at a similar 
fiux for the next 10^ s, before entering a rapid decay 



phase a = 4.60 



+0.18 
-0.14' 



rather similar to GRB 101225A 



and 11 1209 A. This is followed by a shallower decay a = 



0.47 



+0.11 

-0.33' 



before entering a final decay {a = 1.44_Q;Qg) 
after 140 ks, this late time decay is rather typical o f GRB 
X-ray afterglows at these times post-burst (e.g., Evans 
er^[2009). 



In addition to the overall morphology of the X-ray 
emission, a striking feature of the bursts is the rapid 
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apparent dipping and flaring behaviour (some specific 
examples are shown in the inset panels in Figure Isl). At 
these times the emission fades or brightens by a lactor 
> 10 — 100 on timescales of ~ 100 s. These temporal 
changes are unusual in GRB afterglows, but are rather 
typical in GRB prompt emission, which is seen to be 
highly variable, indeed the rapidity of the fiux change, 
At/t ~ 0.01 in several cases, is typical of t he expecta- 
tions of ongoing prompt-like emission (e.g., loka et al. 
|2005j). This change in flux is accompanied by a change 
in the X-ray spectral slope, with brighter periods corre- 
sponding to harder emission, the same general spectral 
evolution a s seen in GRB prompt emission (e.g., Dezalay 
et al.|p^97 ). Finally, if we interpret the early X-ray tlux 



as long-lived prompt GRB emission then the rapid decay 
at the end of this activity could correspond to the obser- 
vation of high latitude emission, as is seen with similar 
decay slopes, but at much earlier ti mes, in many "nor- 
mal" GRBs (e.g., [ Zhang et al.||2006 ). 



2.3. X-ray spectroscopy 

The bright, and long-lived high energy emission pro- 
vides high count rates for the Swift XRT, resulting in the 
possibility of detailed X-ray spectroscopy. However, this 
also represents a challenge in the sense that the early 
X-ray data (< 10^ s) exhibits substantial rapid variabil- 
ity, and any changes in the X-ray spectra associated with 
this will complicate the fitting. We obtained time-sliced 
spectra from the UK Swift Science Data Centre GRB 
Repositor}p| and grouped these such that a minimum 
of 20 counts populated each bin, the time slicing was 
made as fine as possible, following the variations in the 
decay rate, while obtaining a similar number of counts 
per spectrum. Spectra were fitted within Xspec 12.7.1, 
using chi-squared statistics. Galactic column densities 



were taken from the LAB HI Survey ( Kalberla et al. 
2005) and we adopted Solar abundances from "Anders 
Grevesse (1989 1) and cross - sections from Balucinska- 
Church &; McCammon (1992). Intrinsic columns and 
olack body parameters are measured in the source rest- 
frame using the redshifts derived in Section |3] In photon 
counting (PC) mode we fitted all spectra over the en- 
ergy range 0.3-10 keV. In window timing (WT) mode we 
used the energy range 0.4-10 keV for GRB 101225A. In 
the XRT field of GRB 11 1209 A there are a number of X- 
ray sources, and to minimise nearby-source contributions 
to the WT mode spectra we used the energy range 0.5- 
lOkeV. Uncertainties are reported at the 90% confidence 
level. 

In the case of GRB 101225A, the best fitting spec- 
tral model for the early data was claimed to be a black 
body at a temperature of ~ 1 keV, combined with a 



moderate intrinsic X-ray ab sorbing column (Thone et al. 
2011 Campana et al.|[2011 ). However, this model works 



best predominantly at low redshift, and is a somewhat 
worse fit at z = 0.847 (see below). Therefore, it is 
worthwhile investigating the X-ray spectroscopy of both 
GRB 101225A and GRB 111209A in order to ascertain 
what spectral components may be at play, and the con- 
sequences these may have for the nature of the bursts. 
Given the large scale changes apparent in the spectra of 



both bursts we chose to time-slice the resulting data, and 
fit the spectra separately within each time-slice. 

For GRB 101225A the time slices chosen are shown in 
Table [l] At all but the earliest times the spectra are 
adequately fit by a simple absorbed power-law (albeit 
with a modest dispersion in the intrinsic absorption, A/'h, 
likely due to the degree of degeneracy between A/'h and 
photon index), with the resulting fits also shown in Ta- 
ble [l] For the earliest data, which has a fit statistic of 
X^/dof = 193.96/159 for a power-law model, we investi- 
gate additional components. We find that several possi- 
ble models introducing extra parameters can be used to 
improve the resulting fit, this includes the addition of a 
black body component, the use of a cut-off power-law (as 
is often used in prompt GRB-spectra), or the inclusion 
of a break in the power-law, which might be appropriate 
to the early afterglow if a spectral break lies in the X-ray 
band. We discuss the implications of these models below. 

For fits utilising an additional black body, we find 
two broad minima in with temperatures (in the rest 
frame) of 0.08^0.02 and 1.40^0 2 keV respectively (both 
with good xVdof of 152.25/157 and 152.16/157). This 
suggests a significant degree of degeneracy of the in- 
ferred temperature with other parameters, in particular 
the power-law slope and the intrinsic absorption. How- 
ever, we also note that the higher temperature solution 
prefers a relatively low, but poorly constrained intrinsic 
absorption (A^h < 2 x 10^^ cm~^) and hard underlying 
power-law {A{E) ex with, F ^ 1.4 ±0.2), and would 
hence need significant evolution in both to be consistent 
with the later time spectrum. 

A power-law with an exponential high energy cut-off 
also provides an adequate description of the data, which 
is of similar quality to either of the black body fits above 
with x^/dof of 151.92/158, however it also requires a 
low A^H, and very fiat spectral slope with photon index 



F = 0.5^q"^, and a cut-off energy at E^, 



3.1 



+ 1.0 
-0.3 



keV. 



Similarly, a broken power-law fits the break at around 
the same energy as the inferred cut-off above, £^bk = 
3.llQ;5keV, with photon index below and above the 

break of F = 1.20to;l6 and F = 2.1^0.3 respectively, 
again requiring a low TVh- We note that the high energy 
slope of ~ 2.1 is consistent with the inferred slope of the 
early BAT observations F = 1.82 ±0.32, although with a 
large error. The inferred change of slope of AF = 0.9±0.4 
does not provide a strong discrimination of the possi- 
ble nature of the break, a cooling break would imply 
AF = 0.5, which is marginally consistent with the above 
break (although we note that the early time optical spec- 
tral energy dis tributions of these events, in particula r 
GRB 101225A dThone et al.TMTj [Campana et al.|20Tl|) 
show clear deviations from the expectations of a GRB 
fireball model. 

Finally, we also attempt to fit a redshifted ionized 
absorbing model instead of the normal cold absorber. 
While not commonly used for GRBs, such mo dels are 



frequently fit to AGN (e.g., George et al. 1998), and so 



www.swift.ac.uk/xrt_spectra (Evans et al. 2007,2009), pro- 
cessed using the bwift software released in HEASOFT 6.12 



may be appropriate if these events are some class of tidal 
disruption like flare. This also provides a good fit to the 
data, with the resulting fit having x^/dof = 146.46/157. 

We adopt the same approach for our fitting of the 
GRB 11 1209 A XRT spectra, slicing them broadly in 
time to sample both the long-lived plateau, the rapid 
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decay and late-time afterglow-like emission. As with 
GRB 101225A some time slices are inconsistent with 
a simple power-law. Fitting these with an additional 
black-body component provides an adequate fit at tem- 
peratures of either kT = 0.13 or 1.61 keV, with x^/dof 
= 728.75/712 and 703.25/712 respectively for the time 
slice 800-1950 s post burst. Broken and cut-off power- 
law models also provide an adequate fit to the data. 
In the case of a broken power-law the best fit pa- 
rameters (in the same time slice as above) are Fi = 



the second the GRB is c reated by a massive star collapse, 
inside a dense envelope ( Thone et al.|201l" ), at a redshift 
of z ^ 0.33, based on a photometric tit to the late time 
light with a type Ic supernova, and the presence of signif- 
icant X-ray absorption. The interpretation of the origin 
of this burst was comp hcated by its apparently feature 



1.18+^:^^, r2 = 1.59+^;^^, and i?bk = 3.5ll°io keV 
(x^/dof =716.36/712), while for a cut-off power-law we 



+0.03 
-0.03' 



E^^t - 7.65toil keV, with xV^of and^ lThonreff^j^^ et Il.||2QQ9p . 



find F = 0.90 
=681.66/713. 

In summary, for both GRB 101225A and 
GRB111209A, we find that at early times there is 
significant deviation away from a simple power-law 
description of the data. This may imply an unusually 
rapid variation in the underlying power-law spectrum 
that is not fully captured by our time-slicing approach, 
in the sense that the rapid flux variation might require 
the superposition of multiple power-law components, in 
which the impact of the absorption provides a turnover 
at low energies which can be fitted as a black body peak, 
or via other additional components. In this case the 
better fits allowed by more complex models may simply 
be the consequence of adding additional parameters. 
However, we note that cutting down the spectra further 
in time results in fewer counts per spectrum, and so also 
precludes the fitting of more complex models. We note 
that si mply inferr ing the photon index from the hardness 
ratios (Evans et al. 2009) shows a high degree of rapid 
variability, often exhibiting AF ^ 0.5 on timescales of 
only a few tens of secondj^ 

Equally, we cannot rule out the possibility of an ad- 
ditional component in the X-ray spectra, but the data 
do not cleanly distinguish between black-body, cut-off 
power-laws or broken power-laws. We note that in most 
cases the fit statistic for these fits is x^/dof < 1, suggest- 
ing that, if anything, these models over- fit the available 
data, and hence we cannot strongly discriminate between 
different possible physical processes. Equally, since there 
are reasons that each may provide a reasonable physical 
explanation for the observed emission we do not have a 
strong null hypothesis to prefer. 

3. REDSHIFTS 

Critical to any understanding of these events is, of 
course, their distance, something which was missing 
for all GRBs until the discovery of afterglows in the 
late 1990s. While a redshift determination was rela- 
tively straightforward for GRB 11 1209 A, in the case of 
GRB 101225A this has proved particularly problematic. 

3.1. GRB 101225A 

The distance of GRB 101225A has been a matter of sig- 
nificant debate since its initial discovery, with two mu- 
tually exclusive models being proposed to explain the 
event. In the first, the GRB is created by the tidal dis- 
ruption and accretion of an asteroid mass of mate rial 
onto a Galactic neutron star (Campana et al. 2011), in 



less afterglow spectra ([Wiersema et al. 2010 Chornock 
et al.|201Qa Thone et al.|201ip , moderately low Galactic 
latitude (6 = —17''), and a position < 100 kpc in projec- 
tion from M31, if at the same distance as M31. 

However, both pre-imaging, and imaging obtained at 
late times after the burst show that there is a faint, co- 
inci dent quiescent counterpart to GRB 101225A (Fig H 



See http://www.swift.ac.uk/burst_analyser/00509336/ 



obtained Gemini GIVIUS spectroscopy of this source, on 
July 19 & 20 2012. In total a 2.8 hr integration was ob- 
tained in nod and shuffle mode. These spectra were re- 
duced through IRAF in the standard fashion, and clearly 
show two emission lines at wavelengths of 6895 and 9263 
A, as well as further lower significance lines at 9000 
and 9174 A. These lines are naturally explained as be- 
ing due to [Oil], [OIII] and H/3 at a common redshift of 
z = 0.847, thereby finally resolving the problem of the 
distance of GRB 101225A. This is a much larger lumi- 
nosity distance than either of the previously proposed 
models, suggesting that neither fully capture the prop- 
erties of the burst or progenitor. Cut-outs of the 2D 
spectra around the lines are shown in Figure |4j In addi- 
tion to these observations the host of GRB 101225A has 
been observed by the Gran Telescopio Canarias (GTC) 
in July and August 2012. These spectra confirm our red- 
shift of 2; = 0.847, and will be presented in Thone et al. 
(2013, in prep). 

We also re-examined early spectra of the afterglow of 
GRB 101225A with this redshift in hand. These spec- 
tra were obtained from the William Herschel Telescope 
with ACAM on 26 Decembe r 2010 and ISIS on 2 7 De- 
cember 2010, with the MMT ([Schmidt et al.|1989|) on 29 
December 2010, and from Gemini-N with GMOS on 30 
December 2010. A complete log is shown in Table [3j We 
examined each of these spectra individually, and com- 
pared them with a high si^ nal-to-noise ratio compo site 
long- GRB afterglow spectra ( Christensen et al.|2011[ ), as 
shown in Figured While some possible extreniely weak 
features are visible (i.e., as troughs which overlap ab- 
sorption features in the composite afterglow), this com- 
parison suggests that the spectra were not of sufficient 
signal-to-noise for significant features to be seen, unless 
they were particularly strong, and so their non-detection 
is not indicative of a particularly unusual environment. 

3.2. GRB 111209A 

We obtained early spectroscopy of the transient optical 
light of GRB 11 1209 A with both VLT/X-shooter (2011 
Dec 10 at 1:00 UT) and Gemini-N/GMOS-N. These spec- 
tra were processed through the standard pipelines for 
each instrument and show both absorption lines and 
emission lines from the host galaxy, pro viding a red- 
shift oi z = 0.677 ( [Vreeswijk et al. 2011). In addition 
we obtained further spectroscopy from Gemini- S on the 
20 December 2011 and X-shooter on 27 December 2011. 
A log of the ground-based spectroscopic observations of 
GRB 101225A and GRB 111209A is presented in Table[3l 
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and images of the spectra can be seen in Figure jsj 

These data are complemented by two Hubble Space 
Telescope (HST) grism spectra taken on 20 December 
2011 and 13 January 2012, which are described in sec- 
tion [4l 

3.3. GRB 121027A 

We obtained spectroscopy of GRB 121027A with 
Gemini-S/GMOS, and with VLT/X-shooter on 29 and 
30 October 2012. As above fuh details are provided 
in Table [U The GMOS observations were taken under 
poor conditions, and have a limited wavelength range of 
^ 5200 — 7900 A, but do show a strong absorption feature 
at 7770 A, which wa s tentatively interpreted as the Mgll 
doublet at z = 1.77 ( |Tanvir et al ||2012[ ). The X-shooter 
spectra span a much larger wavelength range, and show 
multiple absorption lines from Mg, Fe and Al species at 
a common redshift oi z = 1. 773, confirming our c onclu- 
sion from the GMOS data ( [Kruehler et al.||2012D . We 
note that the X-shooter observations also enable the de- 
tection of emission lines from the host galaxy in the IR 
arm, specifically of [Olll] (4959, 5007), with other fea- 
tures lying in regions of high telluric absorption. A full 
description of these data will be presented in Starling et 
al. (in prep.). 

4. ULTRAVIOLET, OPTICAL AND INFRARED 
PROPERTIES 

4.1. Swift UVOT observations 

In addition to our ground-based observations described 
below, the bright counterparts of both GRB 101225A and 
GRB 11 1209 A were well detected in all of the bands of 
the Swift Ultraviolet and Optical Telescope (UVOT). 

We use the UVOT data for GRB 101225A as given 
in [Thone et~aL] ( |201 1 [ ) . For GRB 111209A, we retrieved 
the data from the U K Science Data Centre. Image mode 
data were processed and analysed using version 3.9 of 
the Swift software, and the nominal 5 arcsecond aper- 
ture. The resulting photometry shows the source to be 
well detected in all bands at early times, and visible in 
white light for ~ 10 days after the initial outburst. The 
UVOT lightcurves are shown in Figure [6l while the de- 
rived photometry can be found in Table]4j 

GRB 121027A was only weakly detected by the UVOT, 
having a white magnitude of 21.55 ± 0.25 in a 353 s 
exposure obtained betwee n 77 s and 1192s post-burst 
(Marshall & Evans 2012). In part this lack of de- 
tections is likely due to the larger luminosity distance 
for GRB 121027A in comparison to GRB 101225A and 
GRB 111209A, while at z = 1.773 the bluest bands are 
also blueward of the Lyman limit. 

4.2. Hubble Space Telescope Observations 
4.2.1. GRB 101225A 

GRB 101225A was observed with HST on 13 January 
2011. At this epoch we obtained a single orbit of ob- 
servation, split between the F606W and F435W filters, 
using the Advanced Camera for Surveys (ACS). A full 
log of the data obtained is shown in Table [5] The data 
were processed through multidrizzle, after correction 
for the effect of pixel-dependent charge-transfer efficiency 
(GTE), and bias striping. In the resulting images we 
clearly detect the counterpart of GRB 101225A with 



magnitudes of F606W=24.60 ± 0.04 and F435W=26.24 
± 0.16. These magnitudes are significantly above the 
host level (see below) in F606W, and somewhat above 
the host in F435W. They imply that the observations 
contain a significant contribution from afterglow light. 

In these data, the source is unresolved, being consis- 
tent with the HST point spread function (PSF) at these 
wavelengths. A subtraction of a PSF reveals no signifi- 
cant residuals at the GRB location, suggesting the host is 
especially compact. We do note that there is a marginal 
elongation visible in the smoothed images (Figure 10), 
but that this is close to the limit of significance, ana is 
also along the direction of the read-axis of the CCD (i.e. 
the direction in which any residual CTE is likely to act). 
At a redshift of z = 0.847, a size less than the HST PSF 
of O.OSarcseconds implies a physical radius of less than 
600 pc. This is a small galaxy, comparable to t he most 
compa ct GRB hosts such GRB 060218 (e.g., ^ Svensson 
et al.||2010aD . 



4.2.2. GRB 111209A 



of 



We obtained two epochs of observations 
GRB 111209A with HST/WFC3. The first was 
obtained on 20 December 2011 and the second on 14 
January 2012. At each epoch we obtained observations 
in two broadband filters, one ultraviolet (F336W) and 
one infrared (F125W). In addition we also obtained 
grism spectroscopy in G102 at each visit, covering the 
spectral range 0.8-1.15 /im. 

The grism spectroscopy was processed through the aXe 
software which was used to background subtract, driz- 
zle and extract spectra of the afterglow at both epochs. 
The resulting grism spectra are shown in Figure u\ (note 
that in the second epoch, the central region has been re- 
moved and interpolated over due to contamination from 
the zeroth order of another star). The spectra show an 
apparent change in the spectral slope, from blue to red 
over the 20 day period between the two epochs of ob- 
servation. This may be due to an underlying supernova 
since the host galaxy light in the IR is apparently well 
below the measured magnitude at the time of our grism 
observations. 

The optical and IR imaging was processed through 
multidrizzle in the standard fashion. At the two epochs 
we measure magnitudes of F336W= 23.76 ± 0.04 and 
25.58 ± 0.15 and F125W = 21.94 ± 0.02 and 22.18 
± 0.15 (all AB-magnitudes). Our ground-based obser- 
vations (see below) imply that at the time of our sec- 
ond HST epoch the ^i-band light was dominated by host 
galaxy, and so we assume the host has F336W(AB) = 
25.78 ± 0.15, whereas we assume the other magnitudes 
are dominated by transient light. 

observations 



4.3. Ground-based optical and in 

4.3.1. GRB 101225A 

We obtained several epochs of optical/nIR imaging of 
GRB 101225A using Gemini-North/GMOS and NIR:(^ 
and WHT/ACAM. These observations were reduced via 
the standard pipelines. Photometric calibration in the g 
and i-bands was perfornied rel ative to the PAndAS sur- 
vey ( McConnachie et al.||2009" ), which covered the same 



Some of these data have been independently reported by 
[Thone et al] poTl] > 
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field. This has the advantage of being taken as a native 
cahbration, and has significantly s maller errors for th e 
comparison stars than reported in Thone et al. (2011), 
especially in the band. For other optica l bands we 
adopt the calibration of Thone et al. f2011) for consis- 
tency, and in the IR calibrate our observations relative 
to 2MASS. An updated photometric calibration for the 
^-band is shown in Table |6] The resulting photometry 
in all bands is shown in Table [71 

The optical (r and z-band) lightcurve of GRB 101225A 
is shown in Figure [6] Beyond one day the z-band can be 
adequately fit with a single power-law with slope of in- 
dex ai = 0.34^0.05 (xVdof = 2.696/6). Fitting a similar 
model to the better sampled r-band light curve provides 
a slope of ar = 0.591q;q2, but the resulting fit is poor 
ix^/dof = 36.92/11), indicative of a de gree of additional 
variability. There is a notable plateau in the r-band be- 
tween ~ 15 — 40 days where the magnitude changed by 
Ar = —0.39 ± 0.19. In the i-band over the same frame 
we see Ai = —0.49 ± 0.09, slightly slower than the decay 
above. Both of these decay rates are very slow for GRB 
afterglows at these late times. A plausible explanation of 
the plateau and apparent reddening is the emergence of a 
supernova, although the data are not strongly diagnostic 
of this (see discussion). We note that the optical after- 
glow appears largely disconnected from the X-ray. There 
is no rapid decay co-incident with the decay observed in 
the X-ray, and the optical continues with a moderately 
shallow slope for many days, long after the source has 
become invisible to the XRT. 

We also obtained a final late epoch of ^-band obser- 
vations with Gemini-N on 1 July 2011. At this time we 
clearly detect a source with g = 26.79 ± 0.14, similar 
to the magni tude inferred from pre-ima^in^ of the field 
in PAndAS (iThone et aL]|2011| ), this is likely the host 
galaxy of GRB 1012237r"^ 

4.3.2. GRB 111209A 

We obtained multiple epochs of imaging in the 
ugrizJHK bands from Gemini North, Gemini South and 
the Very Large Telescope. Gemini images were processed 
as above, while VLT observations were processed through 
the relevant ESO imaging pipeline. The IR observations 
were calibrated against 2MASS stars lying within the 
field of view, while the optical images were calibrated 
against the photometric calibration used by GROND 
(Kann - private communication). For the ii-band we cal- 
ibrate against the F336W iJiST images, and confirm this 
calibration for brighter stars using the UVOT li-band. 
The resulting photometry is shown in Table (Sj and Fig- 
ure [6l 

The bluest bands (namely the UVOT li-band, and 
white light and ground-based ii-band, which have sim- 
ilar central wavelengths) show a moderately steep decay. 
After the first day we find that these bands exhibit a 
slope of Qfwhite = 1.38t^;^^ (xVdof = 9.92/11), rather 
similar to the slope inferred from the X-rays at the same 
time {ax = 1.36^0.05)5 although we note that at early 
times, as with GRB 101225A, there is no strong correla- 
tion between the X-ray and optical light. The similarity 
of the late-time slopes suggests that these bands may 
well be exhibiting standard afterglow behaviour. The 
spectral slope from X-ray to the ii-band at ~ 1 day is 



/^ox ^ 1, while the slope inferred from the X-ray alone is 
Px = 1.39 ±0.07 (Table|2|, these slightly different slopes 
do allow for a cooling break lying between the optical 
and X-ray, although this would be dis-favoured by the 
near identical decay rates. 

In contrast to the bluer bands, the redder optical and 
nIR bands show a markedly slower evolution, resulting 
in a gradual reddening of the afterglow with time. The 
J-band lightcurve, fit with a single power-law from the 
GROND po int at 0.74 days CKann et al.| |2011[ |Kann fc| 
Greiner 2011) until late times results in a best tit de- 
cay of aj = 0.50t^;^4, with a poor xV^of = 39.18/6. 
This suggests both that the decay rate overall is much 
slower than in the i^-band, but also that a single power- 
law is a poor fit to the available data. This poor is 
predominantly caused by the flat behaviour between our 
first i/ST observations and those with HAWK-I ~ 5 days 
later. While colour terms potentially create a problem 
here since F125W and the HAWK-I J-band are not iden- 
tical filters, we do not believe this is significant, since the 
use of secondary calibrators in each field suggests the af- 
terglow was genuinely flat over this period, which in the 
rest-frame corresponds to 7 — 10 days post-burst. The 
difference in the slopes between the J-band and the u- 
band is inconsistent with a simple spectral break between 
the two, since we observe Aa = 0.86 ± 0.07, compared 
to the expectation of the cooling break of Aa = 0.25. 
This is indicative of an additional component contribut- 
ing to the optical/NIR light at late times. In section [t] 
we consider the possibility that this is an accompanying 
supernova. 

To estimate the magnitude of the host galaxy we 
obtained a series of observations from Gemini-S {u^g 
and r-bands on 24 June 2012) and Gemini-N (J-band 
on 6 November 2012). These provide magnitudes of 
u = 25.36 ± 0.25, g = 25.2 ± 0.07, r = 24.94 ± 0.06 
and J> 24.0 (all AB). Given the flatness of the li-band 
between these epochs we assume that all the magnitudes 
measured here are of the host galaxy. However, we note 
that while this seems likely for the i^-band (unless the af- 
terglow has exhibited a prolonged plateau), the g— and 
r-band data remain on the extrapolation of the earlier 
afterglow (see Figure |6|, and so formally these magni- 
tudes provide only an upper limit on the host galaxy 
brightness. 

4.3.3. GRB 121027A 

GRB 101227A was observed with the Anglo Aus- 
tralian Telescope (AAT) beginning 3.5 hr after the burst. 
At this epoch imaging with the Infrared Imaging and 
Spectrograph (IRIS2) revealed the infrared counterpart, 
which was seen to brig hten over the course of the 
first ^ 24 h r post-burst ( [Starling et al. 2012a Levan 



et al. 2012| ). Follow-up observations were obtained 
with Gemini-South/GMOS, showing the source to be 
relatively blue. These observations also demonstrated 
that, as was the case for both GRB 101225A and 
GRB 11 1209 A, the optical and X-ray lightcurves bear 
little resemblance to each other. In particular, over the 
timeframe during which the X-ray lightcurve decays by 
a factor of ~ 10, the optical lightcurve in fact brightens 
by a factor of 3. This comparison and the detailed prop- 
erties of the lightcurves will be discussed in more detail 
by Starling et al. (in prep.). 
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4.4. Longer wavelength observations 

All of the bursts discussed here have been observed 
at millimeter and radio wavelengths by various authors. 
For GRB 101225A, no radio source was detected in ob- 
servations with the Jansky Very Large Array (JVLA) 
befflnning on Dec 30.03 UT to a limit of 60 /iJy at 5 GHz 



(Zauderer et al. 2011a). A second epoch obtained Jan 
6.99 UT also found no source, with limits of 42 /iJy and 
30/iJy at 4.5 and 7.9 GHz respectively (Frail 2011). The 
lack of X-ray detections at this epoch means that only 
optical data is available, and so the resulting broad-band 
SED is poorly constrained. 

Observations of GRB 11 1209 A also initially failed to 
yield any detection, with an upper limit of 132 /iJy at 
34 GHz bein^ placed ap proximately 2 days after the burst 
Hancock et al. 2011). However, observations taken 5 
ays post- burst did reveal a source with flux densities of 
0.85 lb 0.04 mJy (5.5 GHz ), 0.97 ± 0.06 mJy (9 GHz) and 
3.23±0.05 mJy (18 GHz) ( [Hancock et al.|2012|). Our opti- 
cal SED is poorly sampled at arouncTthis time, but aside 
from the steep slope between 9 and 18 GHz, the overall 
shape of the broadband SED is not dissimilar from the 
cano nical exp ectations of a GRB afterglow at this epoch 



(e.g.. Sari et al. 1998). 



Finally; GRB 121027A was observed with the 
APEX/LABOCA bolometer, 4.02 days after the initial 
trigger. These observations failed to yield any sources, t o 
an upper limit of 21 mJy/beam (3(7, |Greiner et al.|2012 ). 



5. HOST PROPERTIES 

The nature of the host galaxies can provide important 
clues to the possible progenitor systems of the ultra-long 
GRB events. For GRB 101225A and GRB 111209A, faint 
underlying host galaxies have been uncovered in deep late 
time imaginj^ 

The constraints on the luminosity of the host of 
GRB 101225A come from late-time imaging with Gemini 
and GTC, and prior im affln^ of the region from the PAn- 
dAS dThone et alTjlMinMcConnachie et al.||2009D . We 
infer that the galaxy is faint and blue, with an absolute 
magnitude Mb ~ —16.3. It is also notably compact, 
being unresolved in ground- and space-based imaging, 
suggesting a size of < 600 pc in radius. The relative sizes 
and absolute magnitudes of these galaxies in comparison 
to other GRB hosts are shown in Figure [Sj 

The host spectrum shows relatively strong, narrow 
emission lines: although the continuum is only very 
weakly detected, we estimate the rest-frame equivalent 
widths EW[oii] ~ 100 A and EW[oiii] ^ 110 A, and find 
them to be consistent with a high specific star-formation 
rate. The lack of absorption lines, particularly of Mgll 
may be suggestive of a relatively small path length of 
cold ISM along the line-of-sight, but the signal-to-noise 
ratio is not sufficiently high for such diagnostics to be 
constraining. 

Given the relatively low signal-to-noise ratio extracting 
measures of metallicity, such as R23 = ([Oil] + [Olll] 
)/ H/3, is challenging, especially since H/3 is in a sky- 
line. However, it is clear that H/3 is not unusually strong 
relative to the strengths of the oxygen lines, and so an 

It is still too early (2.5 months at the time of writing) to 
uncover such a host for GRB 121027A, as it is apparently still 
dominated by the afterglow 



extremely low metallicity is unlikely. 

The host galaxy of GRB 111209A is also faint {Mb ^ 
— 17.6) and compact. Our i/ST imaging at 5 weeks has 
a magnitude consistent with that obtained at much later 
times, suggesting that we are observing the host galaxy, 
but this galaxy is barely resolved in the F336W im- 
age (there is a weak extension visible making the object 
marginally larger than the HST PSF, see Figure 10), 
again suggesting a compact size. Measuring the W% 
light radius from the weak extension seen implies a size of 
~ 700 pc. The VLT/X-shooter spectrum shows emission 
lines from the host. Again, our estimate of the equivalent 
widths of [Oil] and [Olll] are around 30-60 A, although 
this is complicated by the necessity of accounting for the 
contribution of transient light to the latest spectrum. 

Given the better signal-to-noise ratio in this spectrum 
it is possible to estimate i?23, although with some un- 
certainties as to the continuum level. Doing so gives 
log(i?23) ~ 1-0, corresponding to a metallicity of 12 + 
log(0/H) = 8.3 ± 0.3, significantly sub-solar, but not 
especially low. 

The properties of these galaxies, in magnitude, size 
and metallicity are similar to those of compact star 
forming field galaxie s in the same redshift range (e.g. 



Guzman et al. 



1997) and low luminosity blue compact 



dwarf galaxies m the local universe (Figure [9|. In other 
words, they are not unprecedented, but are unusual, and 
suggest a progenitor preferentially found in a dense, low- 
metallicity, intensively star- forming environment. While 
this sample is clearly very small it is interesting to note 
that these galaxies are somewhat unusual, even when 
compared to GRB host galaxies (as seen in Figure [8| - 
they are both smaller, and less luminous than the large 
majority of GRB hosts. Only the host of the low lumi- 
nosity GRB 060218 appears broadly similar, given that 
GRB 0602 18 was also an excep tionally long and atypi- 
cal GRB ( Campana et al.||2006 ) this similarity may not 
be entirely co-incidental. It could suggest a similarity 
in emission mechanisms between the long duration, low 
luminosity bursts, and the ultra-long, but much higher 
luminosity events considered here. 

5.1. Astrometric constraints 

A crucial diagnostic as to the nature of any trans ient 
event is its location within its h ost galaxy (e.g.. B loom] 
eFarp002l 'Fruchter et al."20061 ^enssone^_aLj[MnEf 
Fong et al. 2010; Anderson & James 2009 ^ Events occur- 



rmg at locations mconsistent with the nucleus of the host 
are unlikely to be associated with accretion onto a su- 
permassive black hole, either as AGN outbursts, or tidal 
disruption events (u nless the black hole itself has been 
kicked e.g., |Komossa & Merritt 2008). Nuclear events 
can favour accretion onto supermassive black holes, al- 
though for individual events they do not rule out stellar 
scale events, such as nuclear super novae, wh ich may be 



much m ore common than tidal disruptions (jStrubbe 
Quataert 2009). 



To localise the GRB positions precisely compared to 
their hosts, we compare the astrometry from early ob- 
servations when the afterglow is bright, against the im- 
ages taken at late times when the host dominates. Using 
10 and 8 point sources in common for GRB 101225A 
(using two Gemini frames, taken in December 2010 and 
July 2011) ) and 111209A (using the two HST F336W 
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frames) respectively, we find that the offsets from the 
nucleus of the host galaxies in each case are (0.016 ± 
0.020)'' and (0.011 ± 0.038)'' respectively. In other 
words, both sources are consistent with the nuclei of their 
hosts, and likely lie within 150 and 250 pc of the nucleus 
in each case, as is shown graphically in Figure [lO] How- 
ever, these hosts are also extremely compact, with full- 
width-half-maxima (FWHM) which are, at best, barely 
resolved with HST. Therefore, a significant fraction of the 
total stellar light of these galaxies clearly also lies within 
the error radii associated with the astrometric transfor- 
mations. The positions do not rule out events associated 
with the supermassive black hole, nor do they strongly 
disfavour stellar scale events, more akin to classical long- 
GRBs, which lie at locations consis tent with the nuclei 
of their hosts ^ 10 % of the time ( [Bloom et al]|2002 
Fruchter et"aL][2QQ6l ) . 

6. ENERGETICS 

Based on the redshifts and light curves we can estimate 
the total energy of the bursts. Doing so purely from the 
BAT fiuence is impractical here, due to gaps in cover- 
age of the orbit, and the long lived emission. Therefore 
we integrate the energy based on the c ombined BAT- 
XRT lightcurves (using the method of [ O'Brien et al.| 
|2006), and assuming a power-law interpolation over the 
orbit gaps. Integrating over the first 10^ s (observer 
frame) for each burst yields isotropic energy releases of 
Eiso = 1.20 X 10^2 erg (101225A), Eiso = 5.21 x 10^^ erg 
(111209A) and^;^^^ = 7.00x10^^ erg (121027A). We note 
that given the rapid decay and poor sampling later, this 
number is comparable to the total energy release. The 
corrections to more common bands (e.g., 15-150 keV) are 
dependent on the spectral shape, but we note th at in the 
case of GRB 101225A the BAT spectral slope ([Palmer] 
et al. 2010b) is broadly consistent with that measurea 
by the XK'l (F « 1.8), and suggests a modest correc- 
tion. For GRB 11 1209 A the harder spectral slope results 
in a correction factor making the 15-150 keV energy a 
factor ^ 4 larger, while for GRB 121027A, the slopes 
are softer {Tbat = 1.82 ± 0.09), and as in the case of 
GRB 101225A require a rather smaller correction. 

These energies are relatively common for Swift GRBs, 
lying somewhat below the most extreme Ej ^ ^p cases that 



et al.|2011[[van Eerten fc MacFadyen|2012 l). 
lor canonical iet parameters in simple mo d( 



Nonetheless, 



have been observed 



Racusin et al. 



2008, 2009 



Bloom et al.|2009||Tanvir et al.l2010| ). We notetEat wEiIe 
the late tinie mult i wavelength light curves are poorly 
sampled, and perhaps contaminated by additional com- 
ponents. If we assume that the underlying mechanism for 
the late time emission is similar to that for GRB after- 
glows, then we see no compelling evidence for jet-breaks 
within them (the breaks to a steep decay at ^ 10^ s 
are far too steep for a jet-break). The constraints for 
GRB 101225A are weak, since it was not detected after 
the rapid decay. For GRB 111209A, and GRB 121027A 
the XRT observations continue until ~ 2 x 10^ s, and 
suggest any jet break occurs after this time. In practice 
the data are only weakly sensi tive to breaks c lose to the 
end of these observations (e.g., ^Racusi n et a'L|[2 009), due 
to the limited lever-arm provided at late times. Coupled 
with this, more realistic simulations of GRB-jet evolu- 
tion suggest that for many viewing angles (those not di- 
rectly off-axis), the jet-break is effectively smeare d out, 
making its detection extremely challenging ( ^van Eerten 



jet parameters in simple models (e.g.. Frail 
et al.|200T) it is possible to place limits on the jet openmg 



angle ol Uj > 12 degrees (111209A) and Oj > 10 degrees 
(121027A). The resulting beaming corrected energies are 
E^ > 1.2 X 10^^ erg (111209A) and E^ > 1.0 x 10^^ erg 
(121027), again entirely in keeping wit h the energy out- 
put of man y long duration GRB s (e.g 'Frail et al.ll2 00ll 

. 2oogi r 



Berger et al. 2003; Bloom et al. 2003; Racusin et al 



7. DISCUSSION AND INTERPRETATION 

We are now placed to address the nature of these 
ultra-long GRBs, and from hereon will concentrate on 
GRB 101225A and 111209A, as they are better moni- 
tored than GRB 121027A, and at a redshift at which 
much stronger constraints can be placed on their pro- 
genitors. The nature of GRB 101225A has already gen- 
erated a good deal of controversy. Until this work, it has 
not been cl ear if it belonged t o the extragalactic GRB 
population ( Thone et al.||2011 ), or some new population 
of transients within the Milky Way, perhaps originat- 
ing from ti dal disruptions of asteroids by a neutron star 
( [Ca mpanaT et al. 2011), interestingly one of the models 
tha t was first p roposed for GRB production ( Newma n fc| 
|CoX|,1980^ ^Colgate fc Petsche^l981, ). Below, we discuss 
various possibilities lor the origin ol these extremely long 
GRBs. 

7.1. A core collapse origin? 

It is now clear that the vast majority of long-duration 
GRBs are due to stellar core collapse which also produces 
an accomp anying s upernov a explosion ( see Hjorth et al. 
2003; Sta nek et al.. 2003; H jorth fc Bloom 2012, but see 
also |Fynbo et al.||2Q06t |Gal-Yam et al.||2006f iDella Valle 



et al.||2QQ 6 

examples 



)t al.||2Q06MGal-Y^^ et al.[|2006{ |DeIIa 
Ofek et al.| |2007[ for rare possible co 



counter- 

In each case that has been spectroscopically 



confirmed, the supernova appears as a broad-lined SN 
Ic, suggestive of a stripped and compact progenitor. In 
the standard collapsar model, the central engine must be 
active for sufficient time for the nascent jet t o break-ou t 
of the progenitor (Woosley 1993; Bromberg et al. 12012 ), 
with the 7— ray duration being the lifetime of the jet af- 
ter it has penetrated the stellar envelope. The signature 
of this - a fiat distribution of durations when the 7— ray 
duration is less than the breakout ti me - has been seen 
in Swift hmst^ (Bromberg et al.'2012), directly implying 
that most arise from compact progenitors. This lifetime 
is also natural for material accreting onto the nascent 
black-hole from immediately outside the innermost stable 
orbit. However, these durations are incompatible with 
red supergiant progenitors, which pre-explosion imaging 
ties to the SN IIP (Smartt et al. 2009), since their typical 
radii are 100 to almost 1000 RQ^Any moderately rela- 
tivistic jet would require R^/c> 500 s to tunnel through 
a red supergiant, and so would need a very fine tuned en- 
gine to subsequently create GRBs with durations of only 
a few seconds. Even if accretion driven engines were 
present in a significant fraction of SNe, they would be 
choked by the progenitor envelope, and either no, or a 
rather weak GRB might be seen. Such a scenario has 
been suggested as a plausible explanatio n for the low- 
luminosity GRBs ( Bromberg et al.][2011 ), which, while 
weak in 7— rays, often host SNe very similar t o those 
seen in the prototypical SN/GRB 030329 (e.g., [Hjorth 
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et al.||2003 |Pian et ar| | 2006| ), while populatio n III stars 
mi^ht~aIsQ s utter trom similar constraints (e.g., Nakauchi 



et al. 2012). Therefore, a requirement for jetted events 
from stars of large radius is significant late time engine 
activity that can power the jet. This is likely to arise 
from ongoing (fallback) accretion onto the nascent com- 
pact object. 

At first sight it would appear that these two bursts 
are promising candidates for engines in stars with larger 
radii. This might naturally explain the long durations 
since these larger stars also have a reservoir of mate- 
rial a larger radii from the black hole, which co uld power 
longer lived emission than from a compac t star ( Quataert 
fc Kasen| [20121 [Woosley fc Heger||2012] ). Indeed, such 
an origm has been suggested rece ntly for GRB 11 1209 A 
[Gendre et al.||2012|). The long duration of 



the^vent comfortably provides sufficient time for the jet 
to penetrate a more extended star. In this model the 
rapid variability is not uncommon in GRB prompt emis- 
sion, and the steep decay can be naturally interp reted as 
emission from high latitudes ( Zhang et al.||2006 ). 

Equally, it should be noted, that while larger stars may 
provide a natural explanation of the long durations of 
these event (and are probably more astrophysically com- 
mon), the duration of the burst is intrinsical ly thought 
to be linked to the duration of the engine (e.g., Bromberg 
et al.||20l"2 '), not explicitly the radius of the star. Hence, 



while a shorter duration burst effectively rules out ex- 
tended progenitors, the opposite is not true, these bursts 
could arise from compact stars in which the engine has 
been active for an unusually long period. 

7.1.1. Optical/IR constraints on supernova emission 

At redshifts of z = 0.847 and z = 0.677 for 
GRB 101225A and 111209A respectively, typical SNe 
associated with GRBs will peak around the z- and J- 
bands and will show strong suppression shortward of 
the F-band (3000 A in the rest-frame) due to UV metal 
line blanketing. The light-curves of GRB 101225A and 
GRB 11 1209 A are shown in Figure [6) along with expec- 
tations for a SN 1998bw-like SN. The evolving spectral 
energy distributions of GRB 101225A and GRB 111209A 
are shown in Figure [Tj In both events there is some evi- 
dence of ffattening or re-brightening at late times relative 
to their host galaxy levels, as well as a gradual redden- 
ing in their spectral energy distributions. This may be 
indicative of underlying SNe, since they can ffatten, or 
reverse the afterglow decay, and are typically much red- 
der than the afterglows (which should also decay in a 
broadly achromatic manner). While both the redden- 
ing and ffattening of the light curves of these events is 
clear, it should be noted that this behaviour is not well 
matched in either time-scale or peak luminosity with the 
expectations of SN Ib/c seen in most GRBs. 

If these events are due to core-collapse supernovae, but 
not SNe Ib/c similar to those seen in most other GRBs, 
then the timing of the SNe peak, its magnitude, and the 
underlying spectra become far less well constrained. SNe 
II-P (the progeny of large radii supergiants, Smartt et al. 
[2Q09[ ) can reach peak extremely quickly, remam at a sim- 
ilar magnitude for ~ 100 days, and ex hibit blue spectra 



SNe II-L ap pear to be more akin to SNe Ib/c in l ightcurve 
shape (e.g., Filippenko|1997 Arcavi et al.|2012 ), but also 
show broad tia emission which may leave a detectable 
photometric signature in the band containing Ha. Fi- 
nally SNe Iln, characterized by narrow emission lines 
from circumstellar interactio n can also show extremely 
bright peak magnitudes (e.g.. Smith et al. 2007). In other 
words, the behaviour of the counterparts mignt be very 
different from those typically seen in GRBs, if the under- 
lying supernova is of a different type. 

Distinguishing these possibilities within our relatively 
sparse data is clearly difficult. The late-time spectral 
energy distributions are poorly sampled, with modest 
photometric errors and poorly quantiffed contributions 
from ongoing afterglow emission, and the underlying host 
galaxy (both as a source of additional flux, and of dust at- 
tenuation). Therefore we do not attempt detailed fftting 
of supernova templates, but consider the broad proper- 
ties that would be expected for different SN types. 

For GRB 101225A, the late time SED is clearly incon- 
sistent with being dominated purely by a SN 1998bw-like 
SN at a similar epoch (Figurem), although the peak mag- 
nitude if the late time z-band light was from an SN would 
be rather similar to SN 1998bw. In particular it is far 
too blue at late times, given the redshift of z = 0.847. 
The host galaxy level appears to be well below the blue 
late time ffux, and hence it is likely the result of either 
the associated supernova, or ongoing afterglow emission. 
In the former case the SNe would need to be relatively 
blue, perhaps akin to the early time (~ 15 days) UV 
to optical spectrum of SNe II-P, such as SN 2005cs, or 
2012a^ ^ |Pastorello et al.|2006||Brown et al.|2QQ9| [Bay- 
less et af!]2012p , although it would necessarily be much 
brighter with My ~ —18.5. This is rather brighter than 
the under-lumino us Mr = -15.48 for SN 2005 cs ([PaF 
torello et al .|2006|) or My " ' 

iKA^V ' I . I ' 1 T 



17 for SN 2012aw (Bayless 



et al.| | 201 2 |, although bright S N II-P do exist, e.g., ^IN 



during early times in the plateau (e.g., Filippenko 



1997 



[Arcavi et al. [201 2), with pronoun ced blue to red evo ^ 
liition in the U V throughout (e.g., [Bayless et aL][2012 ). 



1992am ( [Schmidt et al.[[1994[ ). In the case of contin- 
uing afterglow emission, the afterglow would still have 
needed to redden considerably from its earlier behaviour, 
but might do so had the early (extremely blue, possibly 
thermal) emission not been related to the classical after- 
glow emission, but to some additional comp onent such as 
the in teraction with a dense envelope as in lThone et al. 
( [2QTT] ). 

For GRB 11 1209 A, the spectra are clearly different 
from SN 1998bw at similar epochs, both in terms of their 
bulk colours, which appear much bluer, but also in the 
details of our HST grism spectra, which fail to show the 
expected broad features of a SN. These grism spectra do 
show a moderately broad rise around the expected posi- 
tion of the Ha-line. It is possible that this arises from the 
underlying host galaxy, but it could also be the result of a 
H-rich SN. Unfortunately the absence of late time grism 
observations to subtract, and the relatively low S/N of 
the observations preclude detailed study of spectral fea- 
tures, since they are weakly detected, and the host con- 
tribution is unknown. As with GRB 101225A, SN 2005cs 
provides a viable possible model if the late time data are 
dominated by a single supernova component (see Fig- 
ure [t]). Photometrically we note that there is a notable 

35 SN 2005CS and SN 2012aw are amongst few SN II-P with good 
UV coverage via the Swift UVOT. 
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plateau in the light curve of GRB 11 1209 A at approxi- 
mately 10-15 days post burst in the J-band. This would 
be early for the peak of an SN, but if it were associated 
with rising SN emission would imply an unusually lumi- 
nous SN, with My ~ — 21, close to that of the so-called 
superluminous-SNe (e.g G al- Yam 2012) (and only about 
0.5 mag fainter than this at the time of an expected SN 
peak for a SN 1998bw SN). However, we also note that 
the rapid decay from this peak would not be expected in 
these events which tend to evolve more slowly. We also 
note that the same plateau is apparently visible in the 
ix-band, and this may indicate prolonged central engine 
activity is a more likely origin, although then the chro- 
matic behaviour of the lightcurve more generally (i.e., the 
blue to red evolution) could not be naturally explained. 

7.1.2. X-ray constraints on supernova emission 

In principle the X-ray observations might provide a 
strong handle on the nature of supernova emission. Early 
supernovae breakouts can exhibit strong thermal emis- 
sion, which provides both a hallmark of the SN activity, 
and, crucially a radius for the emitting black body. GRB- 
SNe appear to show such signatures, and so the presence 
of a black body component could be a strong indica- 
tor of the burst nature (Starling et al 



l|2012bL 

d bylTh 



Sparre 



tarling 2012). In the model suggested by|'rJione et al 
Q n p , theT)lack body radius inferred is only a few solar 
radii across, and so would appear to immediately rule 
out the collapse of a large radius star. However, this is 
not without difficulties, as we have shown the interpre- 
tation of black body components within the lightcurves 
is far from straightforward, and it is possible that no 
such components exist. Even if we allow a black body, 
the temperature is relatively poorly constrained, with 
both GRB 101225A and GRB 111209A allowing a lower 
temperature black body fit {kT ~ 0.1 keV), as well as a 
higher temperature model. In these low kT cases the in- 
ferred radii are much larger (of order tens of solar radii) , 
and so the constraints allowed on the radius of the pro- 
genitor via X-ray observations are unfortunately weak. 

7.2. A tidal disruption origin? 

Taken at face value the location of both GRB 101225A 
and GRB 11 1209 A consistent with the nuclei of their 
host galaxies suggests a possible link with the super- 
massive black holes which may reside in their cores, 
and could favour a tidal disruption scenari o as was the 
case for GRB lllQ328A/gmft J1644 +57 ("Levan et al. 
^2011; Bl oom et al ||2011| [Burrows et a l. 2011; Zauderer 
et al.,201llD| ), and has also been suggested fo r another ex- 
treme ly lonk transient. Swift 2058.4-h0516 (Cenko et al. 



2012). 

M three bursts (101225A, 111209A and 121027A) 
reach peak luminosities in excess of 10^^ erg s~^, higher 
by an order of magnitude than those seen in Swift 
J1644+57 and Swift J2058+0516. They also have a du- 
ration over which they are detected in 7-rays that is an 
order of magnitude smaller, so they are far from natural 
analogs. In particular, for a main sequence star of mass 
and radius i?*, the circular orbital period (Pt) at the 
ti dal radius (r, ^ R^M bh/M^Y^^) is Pt - I^^M^/Mq 
s (Krolik & Piran 2011). This might be considered the 
minimum timescale on which activity could be observed. 



In practice simulations suggest that the rise time of a 
TDE is at least this, with activity expected for much 
longe r (> 10^ s, e.g., Ayal et al. 2000 Lodato & Rossi 
2011[ ). Although the details of the orbit can create r ather 



ditteren t events, as suggested for Swift J1644+57 (Can- 



nizzo et al. 2011), it would seem challenging to create 
a classical TDE-like event with properties similar to the 
observed bursts. 

However, the high luminosity is significantly in ex- 
cess of the Eddington limit for a 10^^ black hole. 
This, combined with the rapid temporal variations and 
the non-thermal nature of the spectra, suggests t hat any 
emission would be relativisticall y beamed (e.g., Bloom 
et al.|[2011[ [Burrows et"al]|2011| [Zauderer et~aLl[2Ullbp . 



If a jet-component were present then the variations 
could potentially be induced by Lens-Thirring preces- 
sion (Ston e fc Loeb[[2Q1 2), with the rapid cessation at 
~ 10^ s being caused by the precession of the jet out of 
the line of sight. In other words, the rapid apparent end 
to the X-ray emission could be interpreted as an effect of 
jet-precession, rather than of the engine ceasing to func- 
tion (although it by no m eans clear if such je t precession 
can actually occur (e.g., Nixon fc King[[2Q 13). If this 
were the case then the timescale of X-ray activity would 
not rule out main sequence disruptions. In the case of 
GRB 101225A the early optical spectral energy distribu- 
tions are extremely blue, too blue to be accou nted for by 
standard fireball models (e.g., Thone et al.|2011) , in this 
case the optical/UV may be explained by th e presence 
of a hot dis c, as expected in TDE flares (e.g., [Lodato fc 
Rossi 2011[) , which may explain why the optical evolution 
is largely decoupled from the X-ray. In the jetted emis- 
sion scenario, we might expect to observe luminous radio 
emission, as has been the case in other relativistic TDE 



candidates (Zauderer et al. 2011b Cenko et al. 2012) 



The radio upper limits reported in section 4.4 show that 
5-days after outburst GRB 101225A had a 5 GHz flux 
of < 0.06 mJy, at the same time Swift J1644-h57 exhib- 
ited a flux of ~ 2mJy. Even accounting for the differ- 
ence in luminosity distance (a factor of ~ 3) we would 
still expect to easily detect a similar source. In contrast 
GRB 111209A was detected, at a level (~ 1 mJy) which 
is broadly consistent with the extrapolation of the radio 
flux from Swift J1644+57 io z = 0.667. 

If these events are TDE related we might expec t the 
late t i me lightcurves to be similar to the (e. 



1988 



roug. 



Rees 



Lodato fc Rossi[|201 1) expected for TDEs, and 



i ly seen in the case o f Swift J 1 644 -h 5 7 ([Levan et al. 



2011{ iBloom et al.[[2QTT| [Zauderer et al.|[2DT2| . While 
GRB i01225A is not detected after its rapid decay, the 
afterglow of GRB111209A shows a late time slope of 
a = 1.36 ± 0.05, while GRB 121027A shows a slope of 
a = 1.44 ±0.08. These are both close to the t"^/^ slope, 
but are also rather close to the typical late time decay 
rates for GRB afterglows ( Evans et al.[[2 009[), and while 
there is evidence for softening during the rapid decay of 
the lightcurve, the resulting late time spectrum is still 
well fit with a power-law, without the disc (black body) 
emission that might be expected from a tidal event. 

It is interesting to note that the very low luminos- 
ity host galaxies to these bursts could harbour partic- 
ularly low mass black holes (< 10^ M©). In these cases 
the tidal radius of the hole becomes sufficiently small 
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that degenerate stars (in particular white dwarfs) can be 
shredded by the tidal field, rather than being swallowed 
directly. If these events were created by white dwarf dis- 
ruptions then many of the timescale concerns above are 
removed, since the much smaller tidal radius for a dense 
white dwarf leads to an orbital period at disruption of 
Pt ~ 10M^/Mq s, while tightly bound material might 
have an orbital pe riod of several thousand seconds ( Kro- 
lik fc Piran| 2Qll|). This scenario would also naturally 
explain why the host galaxies were extremely low lumi- 
nosity, which is not naturally explained by the main se- 
quence hypothesis. These dense stars might then create 
rather powerful GRBs, that can only occur in low mass 
host galaxies (in contrast to normal tidal fiares which 
shoul d be visible in all but the most massive galaxies 
(e.g., Kesd en|2Q12| . Interestingly such mod els have been 
posit ed before, both for Swift J1644+57 (Krolik & Piran 
20111), long GRBs where no s upernova events are seen 
( |Lu et alT 2008 Gao et al.|2010 ), and previo usly identified 
very long d uration, low luminosity bursts (Shcherbakov 
et al. 2012). However, it remains unclear it such low 
mass galaxies frequently host massive black holes at all. 
The identific ation of Swift J16 44+57 with a LMC-like 
host galaxy (Levan et al. 2011), and the recent identifi- 
cation of a compact X-ray and radio source in Henize 2-10 



long l ived plateau lasting for longer (e.g., [Burrows et al 



2005 Zhang et al. 2006). This naturally means that 



( Reines et al.||201l] ) do imply that some low mass galax- 
ies do harbour such black holes, but their ubiquity (and 
hence the rate at which they may tidally shred stars) 
remains highly uncertain. 

Ultimately, distinguishing between supernova or tidal 
fiare origins is likely to require either an event suffi- 
ciently close that unambiguous supernova signatures can 
be uncovered in its spectrum, or, perhaps more likely, 
the building of a sufficiently large sample of events that 
the locations relative to the nuclei of their hosts can be 
robustly ascertained. 

8. OTHER POSSIBLE EXAMPLES OF THE SAME CLASS 
OF EVENTS 

The previous and existing GRB missions, including 
Swifts are not ideal for detecting very long-duration 
events. Most earlier missions rely on rate triggers, and 
so are at a disadvantage for longer-lived, but lower lu- 
minosity outbursts. Moreover, Swift is in low Earth 
orbit, and therefore most sources are only visible for 
about 45 minutes each orbit. In practice. Swift tends 
to dwell on multiple pointings each orbit, again reduc- 
ing its efficiency for detecting very long events, which, 
of course, typically are detected as "image" triggers. 
However, it is plausible that some very long events can 
trigger the satellite on shorter time scales, for example 
when the lightcurve morphology is such that brighter sec- 
tions of the prompt emission attain sufficient brightness 
to enable shorter image triggers, or rate triggers (e.g., 
GRB 121027A). In any case, the above concerns make it 
clear that relying on a 7-ray duration alone (especially 
from Swift) is not a good route to identifying samples of 
extremely long bursts. 

Given this, we have searched the existing Swift XRT 
repository for examples of X-ray lightcurves which may 
show similar morphology. Such a search is not trivial 
since X-ray lightcurves already exhibit a broad diver- 
sity, and many also show evidence for late time activ- 
ity, typically as ffares found within the ffrst ~ 1000 s, or 



firmly identifying examples will be difficult, although 
there are certain features in the lightcurves of the bursts 
discussed here that do merit attention. In particular, 
three things are striking about the X-ray lightcurves of 
both GRB 101225A and GRB 111209A, and so we have 
used these to ascertain plausible candidates. The cru- 
cial criteria are (i) long lived X-ray emission, lasting to 
~ 10^ s post burst, within two orders of magnitude of the 
peak, (ii) a rapid decay at the end of the plateau, faster 
than the post-jet break expectations, we set a > 3 as a 
constraint and (iii) rapid variations (dips) within the X- 
ray plateau of at least a factor of 5. We therefore identify 
a sample of bursts which meet some of these criteria, and 
based on their lightcurve morphology alone classify them 
as either bronze, silver or gold candidates if they have 1 
(bronze), 2 (silver) or 3 (gold) of the criteria above. We 
note that orbit gaps, and intrinsically faint afterglows 
mean that events which do not show all of the necessary 
features may still lie within a similar morphological class. 

Our search revealed one notable event in the gold 
category, that of GRB051117A (Goad et~al] [20071 ). 
This burst shows bright X-ray emission, m excess of 
10~^ ergs~^ cm~^ for several thousand seconds post 
burst, and also exhibits at least one strong, and nar- 
row dip, lasting for only a few hundred seconds. This 
is followed by a steep break to a decay index of a2 = 
5.3810 40, before settling into a more typical decay of 
0^3 = 0.92 ± 0.04. There is no known redshift for 
GRB051117A, but the resemblance to GRB111209A is 
striking. We also note that GRB 051117A is an extremely 
long burst as measured by the BAT, with a duration 
> 150 s, with some evidence of much later time emission 
dPalmer et al.|2005| . We therefore regard GRB 051117A 
as a likely member of this class of extremely long GRB. 

We also identify several bursts which meet two of 
the criteria outlined above, of particular interest are 
GRBs 111229A, 111016A, 091024 and 060607A. GRB 
060607A lies at z = 3.082 and has tgo = 100±5 s (Ziaeep-| 



our et al. 2008). GRB 091024 has an extreme duration 



of t9n ^ lOUIJs measu red by Swift/B AT and Fermi/GBM 



( Gruber et aL|201 1 [) and was foun d to be at a redshift 
of z = 1.09 fCucchiara et al.||2009|). GRB 111016A has 



oT 

a duration of tgo = 550 ± lOS s, but with possible long 
lived emission beyond the point at which the burst left 
the BAT FOV. GRB 111229A had a shorter prompt du- 
ration, at t go ^ 25 s and was at z = 1.38 (Cucchiara 



et al. 2011|. All of these bursts are apparently more 
distant than either GRB 101225A or GRB 111209A. 

GRB 060607A has marked similarities with 
GRB101225A and GRB111209A, although its gamma- 
ray duration is rather more typical. In the X-ray it 
it exhibits a long-lived plateau persisting roughly an 
order of magnitude fainter than its peak X-ray ffux 
for a duration of ~ 10"^ s, after which it enters a rapid 
decay with a = 3.45^o 12- Iiideed, it would lie in the 
gold category, aside from the lack of obvious dipping or 
ffaring behavior beyond the first few hundred seconds. 
The X-ray and optical afterglow s appear to be largely 
decoupled (Ziaeep our et al.||2QQ8l ), while the host galaxy 
is very faint. Hjorth et al. ~i |2012 ) find the host has 
R > 27.92. The source has also been observed with 
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HST for an exposure of 11194 s in the F775W band with 
ACS/WFC. From this image we find F775W > 30.0. 
This suggests the host galaxy has > —16, this 

is extremely faint, perhaps comparable to the lower 
redshift systems considered here. 

In the case of GRB 091024 the X-ray lightcurve is 
poorly sampled, meaning the overall lightcurve mor- 
phology (and in particular the presence of a steep de- 
cay as in criterion ii) is poorly constrained, however 
there is strong evidence for rapid variability within 
it, which, paired with the very long duration of the 
prompt emission, marks it as a plausible member of this 
class. GRB111016A shows marked flaring, punctuated 
by rapid dips over the first 1000 s, but the Swift orbit 
gap means that the timing (and slope) of a rapid decay 
cannot be ascertained with confidence. In the case of 
GRB 11 1229 A, the plateau is clearly visible, although 
both the level of variability within the plateau and the 
rapid decay at its end (a = 5 ± 3) are more poorly con- 
strained. 

There are several other bursts which meet one of the 
above criteria, for example GRB 0904 17B had a long 
duration, tgo > 260 s, and exhibit ed a bright lat e-time 



flare at about 2000 s post-trigger ( jHolland et al. 2010| ) . 
Thus this population could potentially be moderately 



y 

large, but it seems more likely that the majority of these 
are simply classical long duration GRBs, whose prop- 
erties are outliers to the normal range. This becomes 
more notable when the luminosity of the bursts is con- 
sidered. While GRB051117A lacks a firm redshift, few 
bursts (including those considered above) reach luminosi- 
ties > lO^^ergs"^ several thousand seconds after the 
burst. Since GRB 101225A, 111209A and 121027A do 
attain this luminosity it suggests that if the population 
of ultra-long bursts is significant within the Swift sample 
then the bursts considered here mark the high luminosity 
(and hence potentially most easily detected and studied) 
end of the distribution. 

It is also possible to consider if previous missions have 
identified any members of this class. In particular, Bep- 
poSAX and HETE-2 were able to identify bursts with 
sufficient accuracy for afterglow follow-up, although X- 
ray follow-up typically took place on timescales of sev- 
eral hours at best, such that the X-ray diagnostics de- 
scribed above are of limited value. Nonetheless, the very 
long GRB 020410, detected by B eppoS AX {K\cd.siio et al. 
[2004, ) also exhibited a strong blue to red evolutio n m its 
afterglow a nd arises in a faint underlying host (Levan 
et alj2005D . The HETE-2huTst 021004 exhibits a highly 



unusual early optical afterglow, and arises from a posi- 
tion consistent with t he nucleus of its host galaxy (Fynbo 
let al. 2005; Fruchter et al. 2006), while another burst 
CiRB/XRF 030723A exhibits " a late time bump , and an 
extremely faint host galaxy ( Fynbo et al.||20Q4| . In the 
latter two cases the duration of the pronipt emission is 
apparently much shorter, although the lack of X-ray ob- 
servations limits the strength of any conclusions that may 
be drawn. 

8.1. Rates 

At first sight it would appear that Swift has detected 
at most a handful of these events over its eight year lifes- 
pan, and this would imply that events such as these are 
intrinsically rather rare, a conclusion reached in a sepa- 



rate study by Gendre et al. (2012). However, it is also 
important to consider that observational selection effects 
could act to hinder their detection. In particular, in the 
case of GRB 101225A, the detection was made only by 
integration of the fluence over a long period (over 20 
minutes). Such triggers are typically sensitive only to 
events in which the integrated fluence is rather larger, as 
shown in Figure [ll] The majority of bursts with dura- 
tion > 100 s havenuence > 10~^ergcm~^, roughly the 
median fluence of Swift GRBs. Furthermore, the peak 
fluxes of very long duration bursts are typically lower 
than for the shorter bursts, so that they are less likely to 
result in a rate trigger. In other words there is a clear 
selection bias that Swift cannot detect faint, long lived 
events since they simply fall below the detection thresh- 
olds for any trigger. 

Furthermore, the image trigger durations can be an 
unusually long time for Swift to remain in a single point- 
ing, in particular when long integrations (> 1000 s) are 
needed for detection. To assess this we plot in Figure [12] 
the distribution of Swift snapshot times (Kennea, pri- 
vate communiction) . This shows the total fraction of the 
Swift mission that has been spent in exposures on a sin- 
gle source of > 1088 or 1208 s (i.e. the time needed for 
long image triggers to come into play). The total frac- 
tion of Swift observing time in such exposures in ~ 25%. 
However, even in these cases the trigger is not active 
for the entire exposure time, since an event starting (or 
ending) part way through an observation, but requir- 
ing the full image trigger period to result in a detection 
will not be detected (i.e. if a source that only reached 
the required significance in a 1088s integration were to 
switch on halfway through a 1500s exposure, it would 
not trigger the instrument). Hence, in the extreme, the 
fraction of the total exposure available for the trigger is 
given by the difference between the exposure time and 
the required trigger time (texp — ^trigger)- This number 
can be seen to be rather low (< 10%), although such an 
approach is only approximate, and utilises a rather ex- 
treme example. In practice the true fraction of the total 
Swift exposure time in which it is sensitive to ultra-long 
triggers is likely between 10-25%. For events active over 
several days the limited exposure per orbit is mitigated 
since the exposure can be built up in subs equent orbits 
(e.g.. Swift J2058+0516 |Cenko et a L]|2012|), but for the 
ultra-long bursts (duration of hours) it suggests that the 
true rate to Swift's fluence sensitivity could be a factor 
of several larger. 

The long lived, but typically low peak luminosity 
also suggests that bursts such as GRB 101225A and 
GRB 11 1209 A may be visible over a rather restricted 
horizon in comp arison to class ical GRBs. Indeed, in 



a recent paper Gendre et al. ([2012| suggested that 



GRB 11 1209 A could be detected only out to z ^ 1.4 



significantl y below the mean of GRBs (Jakobsson et al. 



2M2l|2QQ6| ), although the detection of GRB V2m7A at a 
much higher redshift suggests that such events can reach 
higher peak luminosities. Corrected for the limiting vis- 
ibility window and fluence cut it seems likely that the 
rate of ultra-long GRBs to a given fluence may be within 
a factor of a few of that of normal GRBs (with a broad 
range of uncertainty) , and that the observed bias against 
their detection may simply be due to observational se- 
lection effects. If a similar beaming angle were assumed 
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then the true astrophysical rates would equahy imply 
progenitors that were a factor of several rarer than for 
classical long GRBs (although beaming is poorly con- 
strained in these systems). This in itself may provide 
some constraints on progenitor models as they cannot 
be such unusual systems that their rates are significantly 
below those of jet powered coUapsars. 

9. CONCLUSIONS 

We have provided evidence that several ex- 
tremely long-duration transient events, in particular 
GRB101225A and GRB111209A, have very similar 
properties and originated in very similar galaxies. They 
have extremely similar X-ray light curves, luminous 
optical and UV emission and arise from locations close 
to the cores of compact, yet actively star forming 
galaxies. In all of these properties (duration. X-ray 
lightcurve luminosity /morphology, UV luminosity, host 
properties), these GRBs appear as outliers to the 
bulk of the GRB population now being observed by 
Swift. We argue that this likely reflects diversity not 
only in the emission properties of GRBs, but also in 
their progenitors. These bursts have engines which are 
apparently highly active (i.e., luminous 7-ray producing) 
for an order of magnitude longer than is possible in most 
GRBs. This may naturally be explained by the collapse 
of stars of much larger radii than the stripped envelope 
progenitors of most GRBs, which create powerful SN 
Ib/c. Indeed, the late time afterglows of these bursts 
are broadly inconsistent with those of SN Ic, and may 
originate from H-rich supergiants, although evidence in 
this direction remains sparse at present. 

In any case, it is clear that these newly uncovered 
long-duration transients add significant diversity to the 
high energy transient population, either in terms of emis- 
sion mechanisms, or of progenitors (or both). This di- 
verse population may not be intrinsically rare, rather the 
paucity of detection to date may be due to the low sen- 
sitivity of past generations of instruments to long-lived, 
low fiux events. 

Observations to date do not offer a strong confirma- 
tion or rejection of either a core collapse or tidal disrup- 
tion origin for these extremely long GRBs. Clearly local 
events in which SN signatures can be searched for spec- 
troscopically could be extremely important. However, of 
equal importance is the accumulation of a larger sample, 
whose astrometric coincidence with their host nuclei can 
be assessed to determine if they are consistent with an 
origin in supermassive black holes. 
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TABLE 1 

X-RAY SPECTROSCOPIC FITS OF GRB 101225A 
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Note. — The results of time resolved X-ray spec- 
troscopy of GRB 101225A. The timeslices show the time 
since the BAT trigger, in seconds. We show the results 
for simple absorbed power-law fits to the data here, which 
provide an adequate quality of fit (x'^/dof ~ 1) for all but 
the earliest data. With the exception of this bin there is 
no requirement for an additional component to be present 
in the data. Note the A^h is the intrinsic extinction, in 
addition to the Galactic component, and is given in units 
of 10^2 cm-2, z = 0.847. 
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TABLE 2 

X-RAY SPECTROSCOPIC FITS OF GRB 111209A 



Timeslice A^H(int) Photon index x^/dof 



400-800 

800-1950 

1950-2060 

5000-6000 

6000-8000 

10000-20000 

10000-40000 

40000-2.16 X lO^s 



Note. — As for table[T] but for GRB111209A. The table 
shows the results of simple absorbed power-law fits. As can 
be seen, with the exception of the earliest data, these provide 
a generally adequate fit. Note the A^h is the intrinsic extinc- 
tion, in addition to the Galactic component, given in units of 
10^2 cm-2 3t z = 0.677. 



TABLE 3 

Log of ground-based spectroscopic observations of GRB 101225A, GRB 111209A and GRB 

121027A 



Target 


Date-obs 


MJD-obs 


AT (d) 


Telescope 


exp. (s) 


Spectral range (A) 


101225A 


27-Dec-2010 


55557.005 


1.23 


WHT/ISIS 


2400 


5000-9000 


101225A 


29-Dec-2010 


55559.175 


3.40 


MMT/Blue Channel 


4800 


3175-8385 


101225A 


30-Dec-2010 


55560.205 


4.43 


Gemini-N/GMOS 


3600 


3868-6632 


101225A 


19-20 Jul-2012 


56127.506 


571.73 


Gemini-N/GMOS 


10240 


5342-9458 


111209A 


10-Dec 2011 


55905.036 


0.74 


VLT/X-shooter 


4800 


3000-25000 


111209A 


10-Dec 2011 


55905.218 


0.92 


Gemini-N/GMOS 


3600 


3992-8108 


111209A 


20-Dec 2011 


55915.153 


10.85 


Gemini-S/GMOS 


1464 


5992-10000 


111209A 


29-Dec 2011 


55924.116 


19.82 


VLT/X-shooter 


9600 


3000-25000 


121027A 


29-Oct 2012 


56229.253 


1.94 


Gemini-S/GMOS 


2400 


5200-7900 


121027A 


30-Oct 2012 


56230.198 


2.88 


VLT/X-shooter 


6000 


3000-25000 


121027A 


30-Oct 2012 


56230.199 


2.89 


Gemini-S/GMOS 


2400 


5200-7900 



0-39Io.o2 1-20 ± 

0.37lg-°i 1.32 ± 

0.29lg;g^ 1.02 ± 

0.32lg;g^ 1.45 ± 

0.25lg;gi 1.58 ± 

0.29+g-gi 1.68 ± 

0.151^3 1^4^ 

0.19+S;S^ 2.39 ± 



0.01 


700.83/613 


0.01 


908.56/714 


0.02 


408.36/372 


0.02 


385.44/420 


0.01 


526.33/508 


0.01 


604.92/563 


0.05 


103.72/100 


0.07 


62.69/64 



Note. — Ground based spectroscopic observations of our sample of ultra-long GRBs, showing the 
telescopes and instrument setup used for each observation, along with its time since the initial Swift BAT 
trigger. 
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TABLE 4 

Log of UVOT observations of GRB 111209A 



AT (s) 


exp. (s) 


AB-Mag 


Error 


Filter 


742 


19.5 


19.03 


0.26 


UV W 2 


1021 


19.5 


19.87 


0.39 


UVW2 


1369 


19.5 


19.84 


0.38 


UVW2 


5344 


197 


19.66 


0.12 


UVW2 


6780 


197 


19.67 


0.11 


UVW2 


12398 


886 


19.64 


0.06 


UVW2 


19034 


694 


19.93 


0.07 


UVW2 


41800 


886 


20.92 


0.11 


UVW2 


57659 


660 


21.30 


0.17 


UVW2 


75113 


886 


21.75 


0.20 


UVW2 


86682 


886 


21.33 


0.14 


UVW2 


109615 


886 


21.31 


0.15 


UVW2 


123155 


366 


21.56 


0.23 


UVW2 


128450 


886 


21.09 


0.12 


UVW2 



Note. — UVOT observations of 

GRB 111209A, taken in each of the UVOT 
filters. We have attempted to maintain separate 
magnitude measurements in each snapshot where 
possible in order to preserve the time series, but 
note that this occasionally leads to large (> 0.3 
mag) errors on the associated measurements. A 
full log of UVOT observations will appear in the 
online journal. 



TABLE 5 

Log of ifS'T observations of GRB 101225A and GRB 111209A 



Target 



Date-obs 



MJD-obs 



AT 



Filter 



exp. (s) AB-magnitude 



101225A 13-Jan-2011 55574.017 18.241 ACS/F606W 

101225A 13-Jan-2011 55574.034 18.258 ACS/F435W 

111209A 20-Dec-2011 55915.435 11.135 WFC3/UVIS/F336W 

111209A 20-Dec-2011 55915.452 11.152 WFC3/IR/F125W 

111209A 20-Dec-2011 55915.504 11.204 WFC3/IR/G102 

111209A 13-Jan-2012 55939.458 35.157 WFC3/UVIS/F336W 

111209A 13-Jan-2012 55939.455 35.154 WFC3/IR/F125W 

111209A 13-Jan-2012 55939.379 35.079 WFC3/IR/G102 



880 


24.60 ± 0.04 


1020 


26.24 ± 0.16 


1050 


23.76 ± 0.04 


1059 


21.94 ± 0.02 


2212 




1400 


25.58 ± 0.15 


153 


22.18 ± 0.15 


3612 





Note. — A log of the HST observations of each of GRB 101225A and GRB 111209A, showing 
instrument and filters used for each observation, along with the time since burst trigger and where 
appropriate the measured magnitude of the optical/IR counterpart. 



TABLE 6 

B AND 5f-BAND MAGNITUDES OF SECONDARY STANDARD STARS IN 

THE GRB 101225A field. 



REF 


RA 


DEC 


g 


err 


B(AB) 


err 


9 


00:00:48.48 


+44 36:19.3 


19.88 


0.02 


20.89 


0.1 


10 


00:00:47.98 


+44:35:57.8 


19.28 


0.02 


20.12 


0.1 


11 


00:00:50.58 


+44:35:43.5 


19.10 


0.02 


19.88 


0.1 


12 


00:00:51.59 


+44:35:19.1 


19.46 


0.02 


20.02 


0.1 


13 


00:00:43.29 


+44 35:13.1 


18.32 


0.02 


18.89 


0.1 



Note. — Updated magnitudes for ^f-band secondary standard 
stars used in Thone et al. (2011|, and based on the native ^f-band 
calibration from the FAndAS survey. The reference number is the 
ID given in Thone et al. (201 l|. The z -band calibration is broadly 
in keeping with that in I'h one et al.| ([2011), and so we adopt it 
here. Our gf-band calibr ation has sign incantlv smaller errors than 
that reported in Thone ^t al.| ( [2011| , where the calibration was 
extrapolated from yD^y. 
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TABLE 7 

Ground-based photometric observations of the GRB 101225A afterglow 



Date-obs 


MJD-obs 


AT (d) 


Telescope 


band 


exp. (s) 


AB 


-mag 


r 


2(JlU-12-zo 


5555o.947z 


1.1710 


WHT 


r 


300 


22.68 


zt 


n 
u 


08 


zUlU-lz-zo 


55550.9517 


1.1755 


WHT 


i 


300 


23.15 


_l_ 
zc 


n 
u 


25 


zUlU-lz-zo 


5555o.95ol 


1.1799 


WHT 


i 


300 


23.25 


_l_ 


n 
u 


24 


zUiU-iz-ZD 


ooooo.ybUU 


1.8380 


WHT 


r 


300 


22.70 


zt 


n 
u 


07 


2U1U-12-ZO 


5555o.9o43 


1.1880 


WHT 


r 


300 


22.69 


zt 


n 
u 


08 


zUlU-lz-zo 


55550.9083 


1.1921 


WHT 


i 


300 


23.23 


_l_ 
zt 


n 
u 


16 


zUiU-iz-z7 


55557. zool 


1.5100 


Gemini 


r 


120 


22.84 


± 





03 


zUlU-lz-oU 


55560. 18d3 


4.4101 


Gemini 


r 


720 


23.50 


± 





10 


201 1 -01 -OQ 




1 ^ 0^27 


WHT 




1500 


24.10 


± 





14 


2011-01-09 


55570.8500 


15.0732 


WHT 


B 


3600 


25.88 


± 





40 


2011-01-10 


55571.2056 


15.4293 


Gemini 


9 


1200 


25.22 


± 





10 


2011-01-10 


55571.2229 


15.4467 


Gemini 


r 


1440 


24.32 


± 





03 


2011-01-10 


55571.2444 


15.4682 


Gemini 


i 


1890 


24.13 


± 





04 


2011-01-19 


55580.8419 


25.0657 


WHT 


r 


3600 


25.01 


± 





32 


2011-01-23 


55584.2674 


28.4911 


Gemini 


r 


900 


24.87 


± 





08 


2011-02-03 


55595.2632 


39.4870 


Gemini 


9 


900 


26.25 


± 





32 


2011-02-03 


55595.2486 


39.4724 


Gemini 


r 


900 


24.71 


± 





19 


2011-02-03 


55595.2347 


39.4585 


Gemini 


i 


900 


24.62 


± 





08 


2011-02-03 


55595.2181 


39.4418 


Gemini 


z 


1260 


24.90 


± 





25 


2011-07-01 


55743.5576 


187.7810 


Gemini 


9 


3000 


26.79 


± 





14 



Note. — Photometric observations of GRB 10 1225 A obtained at the Wihiam 
Herschel Telescope and Gemini-N. The magnitudes have not been correcte d for 
Galactic extinction, to do so would require corrections based on Schla fly Sz\ 
[Finkbeiner| ( |2011| ) of As=0.37, Ag=0.33, Ar=0.23, Ai=0.17 and A^=0.13. 



TABLE 8 

Ground-based photometric observations of the GRB 111209A afterglow 



Date-obs 


MJD-obs 


AT 


Telescope 


band 


exp. (s) 


mag 


(AB) 


2011- 


-12- 


-09 


55905 


2060 


0.9059 


Gemini-N 


9 


3 X 60 


20 


24 


± 





04 


2011- 


-12- 


-09 


55905 


2025 


0.9024 


Gemini-N 


r 


3 X 60 


19 


82 


± 





02 


2011- 


-12- 


-09 


55905 


1991 


0.8990 


Gemini-N 


i 


3 X 60 


19 


43 


± 





02 


2011- 


-12- 


-09 


55905 


1956 


0.8955 


Gemini-N 


z 


3 X 60 


19 


52 


± 





02 


2011- 


-12- 


-19 


55915 


0757 


10.7756 


Gemini-S 


u 


8 xl80 


23 


41 


± 





08 


2011- 


-12- 


-25 


55920 


0293 


15.7292 


VLT/FORS2 


u 


4 X 300 


23 


97 


± 





10 


2011- 


-12- 


-25 


55920 


0462 


15.7461 


VLT/FORS2 


i 


8 X 120 


22 


36 


± 





03 


2011- 


-12- 


-25 


55920 


0610 


15.7609 


VLT/FORS2 


z 


8 X 120 


22 


18 


± 





08 


2011- 


-12- 


-25 


55920 


0252 


15.7251 


VLT/HAWKI 


J 


35 X 60 


20 


70 


± 





04 


2011- 


-12- 


-29 


55924 


0348 


19.7347 


VLT/FORS2 


9 


4 X 250 


23 


22 


± 





02 


2011- 


-12- 


-29 


55924 


0484 


19.7483 


VLT/FORS2 


R 


4 X 250 


22 


86 


± 





02 


2012- 


-01- 


-02 


55928 


0362 


23.7361 


VLT/FORS2 


u 


6 X 300 


24 


50 


± 





23 


2012- 


-01- 


-02 


55928 


0264 


23.7263 


VLT/HAWKI 


J 


35 X 60 


20 


99 


± 





04 


2012- 


-01- 


-09 


55935 


1356 


30.8355 


VLT/HAWKI 


J 


35 X 60 


21 


32 


± 





10 


2012- 


-01- 


-14 


55940 


0406 


35.7405 


VLT/FORS2 


u 


4 X 300 


25 


29 


± 





24 


2012- 


-01- 


-14 


55940 


0173 


35.7172 


VLT/HAWKI 


J 


35 X 60 


21 


22 


± 





05 


2012- 


-01- 


-19 


55945 


0319 


40.7318 


VLT/FORS2 


u 


9 X 300 


25 


22 


± 





24 


2012- 


-01- 


-19 


55945 


0241 


40.7240 


VLT/HAWKI 


J 


35 X 60 


21 


54 


± 





06 


2012- 


-01- 


-21 


55947 


0652 


42.7651 


VLT/FORS2 


9 


4 X 250 


24 


16 


± 





04 


2012- 


-01- 


-21 


55947 


0787 


42.7786 


VLT/FORS2 


R 


4 X 250 


23 


62 


± 





03 


2012- 


-01- 


-21 


55947 


0351 


42.7350 


VLT/FORS2 


i 


8 X 120 


23 


05 


± 





05 


2012- 


-01- 


-21 


55947 


0499 


42.7497 


VLT/FORS2 


z 


8 X 120 


22 


84 


± 





11 


2012- 


-01- 


-24 


55950 


0425 


45.7424 


VLT/FORS2 


u 


9 X 300 


25 


44 


± 





20 


2012- 


-01- 


-24 


55950 


0275 


45.7274 


VLT/HAWKI 


J 


35 X 60 


21 


56 


± 





06 


2012- 


-06- 


-24 


56102 


4181 


198.118 


Gemini-S 


u 


5 X 300 


25 


36 


± 





25 


2012- 


-06- 


-24 


56102 


3876 


198.088 


Gemini-S 


9 


5 X 180 


25 


72 


± 





07 


2012- 


-06- 


-24 


56102 


4012 


198.101 


Gemini-S 


r 


5 X 180 


24 


94 


± 





06 


2012- 


-11- 


-06 


56237 


2979 


332.998 


Gemini-N 


J 


45 X 60 


> 24.0 







Note. — Photometric observations of GRB 111209A obtained at Gemini-N, Gemini- 
S and the Very Large Telescope. The magnitudes have not been corrected for Galac- 
tic extinction, to do so would require corrections based on Schlafly &; Finkbeinerl 
(12011 ) of A^,=0.08, Ap=0.06, Ar=0.04, A^=0.03 and A;,=0.02, Aj=0.0013, AH=^m^ 
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Fig. 1. — The spectral-hardness (ratio of fluence in 50-100 keV over 20-50 keV) versus duration diagram for CGRO/BATSE GRBs (red 
points) and Swift GRBs (blue points), with the locations of ORB 101225A, GRB 111209A and ORB 121027A marked (note these are 
approximate due to the lack of Swift orbit coverage). These three events have durations much longer than any seen by BATSE. In the 
case of GRB 101225A, the long-lived, low level emission could easily have been missed, while GRB 111209A was seen as an extremely long 
burst by Konus- Wind. 
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Fig. 2. — Parameter space for transients in the 7-ray sky, showing the duration of the burst, and the approximate average luminosity 
over that duration. At low luminosity there are numerous Galactic sources that we do not include in further detail; at higher luminosity 
the outbursts for soft-gamma repeaters (SGRs) in our own Galaxy are shown, as well as extragalactic transients such as long and short 
duration GRBs (LGRBs and SGRBs), and the likely population of low luminosity GRBs (LLGRBs). Two recently discovered very long 
transients, thought to be from tidal disruption events are also shown (labelled TDEs?). The bursts considered in this paper (GRB 101225A, 
GRB 111209A and GRB 121027A) are clearly outliers to any of these aforementioned classes. 
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log ^rest [s] 

Fig. 3. The rest-frame X-ray lightcurves of GRB 101225A (red), GRB 111209A (green) and GRB 121027A (blue), in luminosity space 
(rest frame 0.3-10 keV). There are compared with two well known very long GRBs (namely 060218 (orange) and 100316D (grey), and all 
Swift bursts with known redshift (grey background shading). In light blue we also show the light curve of Swift J1644+57 (Levan et al. 2011 
[Burrows et al. 2011), which persists for even longer, but is rather less luminous. The light curves of the ultra-long events are extremely 
similar, both in overall shape, and the strong dipping behaviour seen. The upper panels show the strong substructure present in the light 
curves of the ultra-long bursts, demonstrating rapid dipping, with relatively narrow dips (in particular in the case of GRB 111209A). We 
also note that the luminosity of these events at late times (~ 10^ s) is substantially (a factor of ~ 100) greater than that of typical GRBs. 
This is likely because the prompt emission in these cases persists for much longer, but clearly marks these events apart from other GRBs. 
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Fig. 4. — Observations of GRB 101225A. [Upper left panel] shows a finding chart of the region around the host galaxy (within blue square). 
The image used is the late-time ^f-band observation made with Gemini-N/GMOS (north up and east left). [Upper right panels] shows our 
best afterglow spectroscopy, overlayed with a high-S/N summed GRB afterglow template (red line from Christensen et al. (20li|). The 
top sub-panel is the WHT/ISIS data (the gap is between the blue and red arms) and the bottom one the Gemini/ GMUy data, in neither 
case are absorption lines significantly detected, but this is not surprising given the signal-to-noise ratio. [Lower left and right panels] show 
cut-outs of our late time Gemini/GMOS 2D spectroscopy of the host galaxy around the prominent oxygen emission lines (circled). 
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Fig. 5. — Finding chart and spectroscopy for GRB 111209A. The counterpart is shown in the blue box on a HST F125W image. The 
right spectrum shows a portion of the absorption spectrum obtained by X-shooter 0.74 days after the burst, demonstrating strong hues of 
Mg and Fe in the ISM of the host galaxy. The lower panels show the 2D spectrum from X-shooter on 29 December 2011, and show strong 
emission lines from Oil and Olll (note that there is still a contribution from the afterglow, as seen in the continuum in these 2D spectra). 
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Fig. 6. — The UV and optical light curve of GRB 101225A (top) and GRB 111209A (bottom). The colours represent the same filters 
in each panel as shown by the inset key. We note that in the GRB 11 1209 A panel we have represented the UVOT white light filter with 
the same colour (black) as the w-band given the similar central wavelength, but where the li-band is indicated with circles, the white light 
points are marked as triangles. In addition we show the inferred temporal slopes as solid lines for the r and i-bands for GRB 101225A, and 
for the u- and J-bands in GRB 111209A. The curves show the expectation of a SN 1998bw-like SN in the relevant band (colours coded as 
above) with no scaling or stretching. This is seen to be similar to the late time magnitude of GRB 101225A, but well below the level seen 
in GRB 111209A 
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Fig. 7. — Evolution of the broad-band spectral energy distribution (SED) of the afterglow of GRB 101225A (top) and GRB 111209A 
(bottom) over several epochs, as indicated. Spectroscopic observations of GRB 101225A with Gemini-N (grey) and of GRB 11 1209 A with 
the HST WFC3 grism are also plotted (at epochs 11 and 35 days post bursts, colour coded for the SED at that epoch). Both events show 
a blue to red evolution, although this is particularly strong in the case of GRB 101225A, which at early times shows a strongly rising SED 
towards the blue, possibly indicative of a thermal component peaking in the EUV bands (Thone et al. 2011). At late times an SN 1998bw 
template (magenta, chosen to approximately match the same rest frame epoch as the GRB data) does not provide a good fit to the global 
SED of either burst, suggesting that the supernova underlying each could be rather different (SN 2005cs is shown in black as a template 
which may represent the data), alternatively there could be a still important contribution from any afterglow emission. 
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Fig. 8. — Constraints on the supernovae and host galaxies associated with GRB 101225A and GRB 111209A, based on the observations 
with the CFHT, Gemini and H ST. The top panel shows the absolute magnitudes of the host galaxies (black, Svensson et al. 2010b) and 
supernovae (red, |Hjorth Bloom 2012) associated with GRBs. The filled triangles represent firm spectroscopic associations with yiNe, and 
the open circles weaker spectroscopic, or photometric SNe. The redshifts of GRB 101225A and GRB 11 1209 A are indicated by dashed 
blue lines, while the open red boxes indicate the approximate magnitudes of the associated supernovae, if the plateau in each lightcurve is 
ascribed entirely to supernova emission. The lower panel shows the radii of several GRB hosts, compared to the limits on GRB 101225A 
and GRB 111209A based on the non-detection (or in the case of GRB 111209A marginal detection) of extension in the //aS'T observations. 
As can be seen, the hosts of GRB 101225A and GRB 111209A are extre me in comparison with normal GRB hosts. However similarly 
compact galaxies, such as IZwl8 ( ^Fiorentino et al.^20ld] [Aloisi et al.|1999| > can be found in the local Universe, and so these properties are 
not unprecedented. 
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Fig. 9. — The properties of the two host galaxies, absolute blue magnitude versus equivalent- width of [Oil] (GRB 101225A in blue and 
11 1209 A in red), in comparison to a sample of local blue compact dwarf galaxies from Kong et al. ("2002) (black triangles), moderate redshift 
galaxies in the GOODS field (open circles, from Kobulnicky Sz Kewley 2004), and a sample of pre- tlwift GRB host galaxies (green crosses). 
Although the equivalent- width measures, in particular, have rather high uncertainty, they are consistent with being members of the blue 
compact dwarfs, but are somewhat offset from the field galaxy sample of Kobulnicky & Kewley (2004), or the GRB host galaxies (although 
all sam ples may contain significant selection effects. The GRB host galaxy Mb values have been tit using the techniq ues of Perley et al.] 
(|2013|>/and the [Oil] equivalen t widths are taken from Bloom et al. (2001 1998); Diorgovski et al. (2001 1998); Le F loc^h et al.| ( |2002| >; 
|(Jasto)-Tirado et al., ( ,2001j ; .Piro et al.| ( [2002} ; [Price eTal., ( ^2UU2, ,2UU3j ; ,Margutti et al., ( ,2UU7j ; ,Frochaska et al., ( ,2004^ ). 
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Fig. 10. — Astrometric constraints on the locations of GRB 101225A and GRB 111209A. In each case the circle shows the location, and 
associated la error on the position of the afterglow, overlayed on a late time image of the host galaxy of the GRB. The cross marks the 
optically determined centroid of the galaxy, determined by fitting a gaussian profile to the host light (in each case the hosts are at best 
barely resolved by the instruments (GMOS for GRB 101225A and WFC3 for GRB 111209A), and so host morphology is not important to 
the determination of the centroid. 
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Fig. 11.— Constraints on th e detectability of ultra-long GRBs with Swift. The lower panel shows the duration fluence relation (see also 
|Levan|2012[|Gen dre et al.'2012), along with the median values in different duration bins (red dots), clearly showing that the longer duration 
bursts have on average a greater fluence, the lack of low fluence longer events suggests that this is an effect intrinsic to the detector (i.e 
there is a selection against faint, long lived transients). The upper panel shows the relation between peak-flux and duration (again red 
dots indicate median values), this suggests that on average longer bursts have fainter peak fluxes, and so are less likely to result in a rate 
trigger for an instrument such as the Swift BAT. 
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Fig. 12. — The distribution of Swift snapshot duration (i.e., the duration between slews) from 2012, the sohd black line shows the 
cumulative fraction of the Swift lifetime spent in exposures longer than the snapshot duration, while the red and blue lines show the 
fraction of time where the longest image triggers (1088 and 1208 s) can be used, which is a very small fraction of the overall mission. This 
suggests that while ultra-long bursts are observationally rare, this may be largely due to a strong selection function. 



